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ABSTRACT 


A  hyperbaric  chamber  environment  is  simulated  using  a 
FORTRAN  computer  program.  This  program  allows  the  user  to 
input  data  files  of  varying  chamber  conditions  and  provides  a 
detailed  output  which  can  be  used  for  design  or  chamber 
operations.  Specifically,  this  program  determines  life 
support,  temperature  and  humidity  control,  and  power 
requirements.  After  developing  the  theory  and  describing  the 
program  in  detail,  complete  examples  simulating  a  hyperbaric 
chamber  at  2000  fsw,  a  submarine  lockout  chamber  at  400  fsw, 
and  a  shallow  underwater  habitac  at  50  fsw  are  demonstrated. 
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CHAPTER  I 


INTRODUCTION 


A  hyperbaric  environment  exists  in  a  facility  which  may 
be  pressurized  from  one  to  several  hundred  atmospheres. 
Edmonds  (1976)  defines  a  hyperbaric  chamber  as  an 
environmental  pressure  vessel  where  the  ambient  pressure  can 
be  increased  from  one  atmosphere  to  many  times  atmospheric 
pressure  depending  on  the  structural  capabilities  of  the 
vessel  concerned.  Hyperbaric  chambers  are  also  known  as 
recompression  chambers  and  decompression  chambers.  A 
recompression  chamber  is  used  for  compressing  a  patient  or 
diver  as  part  of  a  therapeutic  regime  or  for  training.  A 
decompression  chamber  is  used  to  decompress  a  subject  already 
exposed  to  increased  pressure.  Miller  and  Koblick  (1984) 
define  an  underwater  habitat  as  a  facility  which  acts  as  a 
survival  shelter  allowing  humans  to  survive  and  live  on  the 
ocean  floor  continuously  for  several  weeks.  Hyperbaric 
chambers,  underwater  habitats,  diving  bells,  and  submarine 
lock-out  chambers  require  unique  design  considerations  to 
account  for  the  additional  pressures  and  the  effect  on  human 
physiology  and  thermodynamic  properties. 

History 

Hyperbaric  chambers  have  had  medical  uses  for  over  300 
years.  Edmonds  (1976),  Miller  and  Koblick  (1984),  Shilling  et 
al.  (1976),  and  U.  S.  Navy  Diving  Gas  Manual  (1991)  provide  a 
chronological  history  of  hyperbaric  chambers; 
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1662  “  The  first  use  of  a  hyperbaric  chamber  in  medicine  was 
made  by  a  British  physician.  Dr.  Henshaw.  This  chamber  was 
pressurized  and  depressurized  by  a  pair  of  organ  bellows. 

1830's  -  Three  Frenchmen,  Junod  (1834),  Tabarie  (1838),  and 
Pravaz  (1837)  renewed  interest  in  hyperbaric  chambers. 
Junod 's  chamber  was  a  1.6  ro  diameter  sphere  capable  of 
pressures  up  to  4  atm.  He  noted  the  euphoric  effect  of 
compressed  air.  Tabarie  claimed  improvement  in  the  treatment 
of  lung  disease.  Pravaz  built  a  large  chamber  capable  of 
holding  12  patients  at  a  time. 

1850 's  -  A  rapid  rise  in  the  construction  of  hyperbaric 
chambers  was  observed.  Compressed  air  was  provided  by  steam 
engines  and  chambers  capable  of  simultaneously  treating  up  to 
50  patients  were  built. 

1860  -  The  first  hyperbaric  chamber  in  North  America  was 
built  in  Ashawa,  Canada. 

1891  -  The  first  chamber  in  the  United  States  was  built  in 
Rochester,  NY. 

1921  -  Orval  J.  Cunningham  of  Kansas  City,  MO  constructed 
the  world's  largest  hyperbaric  chamber.  It  was  5  stories 
high,  had  a  diameter  of  20  meters,  and  was  equipped  with 
amenities  suitable  to  a  first  class  hotel. 

Since  1920 's  -  Numerous  hyperbaric  chambers  have  been 
designed  for  treatment  of  divers.  with  the  progress  of 
hyperbaric  and  diving  medicine,  chambers  are  becoming 
increasingly  complex. 

1960 's  -  The  French  and  Americans  began  work  on  underwater 
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habitats.  The  French  built  Conshelf  I  &  II  while  the  United 
States  built  Man-in-the-Sea  I  &  II  and  Sealab. 

Modern  Hyperbaric  Facilities  include; 

1)  Ocean  Simulation  Facility  (OSF)  -  located  at  the  U.  S. 
Navy  Experimental  Diving  Unit  in  Panama  City,  FL.  This 
chamber  is  one  of  the  world's  largest  man-rated  hyperbaric 
facilities  with  a  depth  capability  of  2250  fsw.  A  wide  range 
of  experiments  in  biomedical  studies,  human  factors 
verification,  and  deep  diving  systems  performance  tests  are 
conducted  in  wet  and  dry  chambers. 

2)  Naval  Medical  Research  Institute  (NMRI)  -  located  in 
Bethesda,  MD.  This  complex  is  the  U.  S.  Navy's  largest 
hyperbaric  complex  fully  devoted  to  research  and  contains 
three  dry  chambers  and  a  vertical  wetpot  rated  to  2280  fsw  and 
two  dry  chambers  rated  to  3370  fsw.  NMRI  conducts  hyperbaric 
research  focused  on  the  biomedical  and  physiological  problems 
encountered  by  the  fleet  diver. 

3)  Duke  University's  Hyperbaric  chamber  -  located  in 
Durham,  NC.  This  facility  is  capable  of  employing  any  desired 
gas  at  high  pressures  ranging  from  150,000  ft  of  altitude  to 
1000  fsw.  Research  is  conducted  on  acoustic  speech 
modification,  high  pressure  nervous  syndrome  and  carbon 
dioxide  effects. 

A  typical  experimental  hyperbaric  facility  is  shown  in  figure 

1.1. 

Several  factors  must  be  considered  when  designing  a 
hyperbaric  chamber.  NAVFAC  DM-39  (1982)  addresses  the 
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following  design  requirements: 

1)  Pressure  vessel  design  including  material  selection, 
wall  thickness,  welding  requirements,  connections,  view  ports, 
and  doors. 

2)  Fire  protection  and  extinguishing  systems. 

3)  Electrical,  utility,  and  communication  systems. 

4)  Life  support  systems  including  mixed  gas  breathing 
systems,  carbon  dioxide  control,  and  oxygen  replenishment. 

5)  Temperature  and  humidity  control. 

6)  Power  requirements. 

The  objective  of  this  program  is  to  determine  life  support, 
temperature  and  humidity  control,  and  power  requirements. 

During  most  hyperbaric  chamber  designs,  these  factors  are 
analyzed  ^.eparately  and  compiled  into  a  finished  design.  This 
program  allows  the  user  to  input  varying  design  considerations 
and  provides  a  detailed  output  of  design  requirements.  The 
most  unique  capability  of  this  program  is  the  use  of  a 
computer  to  determine  the  thermodynamic  and  psychrometric 
properties.  Traditionally,  psychrometric  charts  have  been 
used  to  determine  dehumidification,  moisture  removal,  and 
ultimately  energy  requirements.  Although  specially  developed 
psychrometric  charts  are  available  for  the  higher  pressures 
required  in  a  hyperbaric  environment,  they  are  time  consuming 
and  lack  the  accuracy  of  a  computer  program.  Additionally, 
reading  values  off  these  charts  requires  non-linear 
interpolation.  . lerefore,  the  necessary  relationships  were 


developed  to  replace  the  psychrometric  chart  and  simulate  a 


hyperbaric  chamber  environment. 


INNER 

DECOMPRESSION 

CHAMBER 


OUTER 

DECOMPRESSION 

CHAMBER 


2. 


3 


Divers  Swimming  and  Working  in  Wetpot  are 
Supported  by  Tenders  in  Igloo. 

Igloo  IS  Primary  Working,  Living  and  Resting 
Space  "Dry”  Dives  are  also  Conducted  in 
this  Space. 

Decompression  Chamber  is  used  for  Sleeping 
and  for  Treatment  of  Decompression  Sickness 
Each  Space  in  the  Facility  Can  Be  Pressurized 
Individually  or  in  Combination  With  Other 
Spaces 


Figure  1.1 


Typical  Experimental  Hyperbaric  Facility 
(Shilling  et  al.  1976) 


Objectives 

Based  upon  relevant  aspects  of  human  physiology,  this 
program  determines  the  ideal  diving-gas  mixture  along  with  an 
acceptable  rai  c.e  of  diving-gas  mixtures  for  a  given  depth  or 
pressure.  Additionally,  the  required  ventilation  rate  without 
carbon  dioxide  absorbers  or  the  amount  of  absorption  material 
(sodasorb)  for  saturation  diving  is  calculated.  Various 
properties  of  the  selected  diving-gas  mixture  are  determined, 
end  by  using  thermodynamic  properties  of  the  gas,  the  heating 
and  cooling  requirements  are  calculated.  Finally,  the  energy 


a 

requirements  for  dehumidification  and  reheating  to  a  desired 
^  temperature  are  determined  along  with  moisture  temoval 

requirements . 
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THEORY  FOR  DESIGN  OF  HYPERBARIC  ENVIR0NNENT8 

Basic  Phvaiolocrv 

A  basic  knowledge  of  human  physiology  is  required  to 
determine  an  acceptable  diving-gas  mixture  for  varying  depths 
or  pressures.  While  a  given  mixture  may  be  acceptable  at  one 
pressure,  it  can  cause  discomfort  or  even  be  fatal  at  other 
pressures.  The  percentage  of  oxygen  required  to  sustain  life 
decreases  as  pressure  increases  in  order  to  maintain  an 
acceptable  oxygen  partial  pressure. 

The  icspiratory  system  of  the  body  is  responsible  for 
oxygen  uptake  and  carbon  dioxide  removal.  Oxygen  uptake  is 
governed  by  two  main  factors  listed  oy  U.  S.  Navy  Diving  Gas 
Manual  (1971); 

1.  Respiratory  Minute  Volume 

2.  Effectiveness  in  which  the  lung  extracts  oxygen  from 
alveolar  gas. 

Respiratory  minute  volume,  which  is  defined  as  the  product  of 
the  respiratory  rate  (breaths/min)  and  the  tidal  volume 
(volume/ breath) ,  determines  the  duration  of  exposure  of  blood 
to  alveolar  air.  The  effectiveness  in  which  the  lung  extracts 
oxygen  from  alveolar  gas  depends  upon  1)  composition  and 
volume  of  air  supplit-d  and  2)  partial  pressures. 

The  transfer  of  gases  takes  place  due  to  differences  in 
partial  pressures.  The  transfer  of  oxygen  occurs  with  no 
complications  if  the  oxygen  partial  pressure  is  between  0.16 
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and  1.6  atm  for  short  periods  of  time.  However,  if  the  oxygen 
partial  pressure  is  below  0.16  atm,  anoxia  (oxygen  shortage) 
which  may  lead  to  fainting  and  possible  brain  damage  can 
occur.  Likewise,  if  the  oxygen  partial  pressure  exceeds  1.6 
atm,  oxygen  poisoning  which  may  lead  to  convulsions  and 
unconsciousness  may  occur.  NOAA  Diving  Mantial  (1991)  states 
that  a  1.6  atm  oxygen  partial  pressure  can  only  be  tolerated 
for  up  to  45  minutes  and  that  the  maximum  oxygen  partial 
pressure  for  long  periods  should  not  exceed  0.5  atm.  NOAA 
Diving  Manual  (1991)  also  recommends  that  for  long  term 
exposures  in  a  hyperbaric  chamber  or  underwater  habitat  the 
oxygen  partial  pressure  of  the  diving-gas  should  be  between 
0.3  and  0.4  atm.  Figure  2.1  shows  the  percentage  of  oxygen  in 
the  diving-gas  mixture  as  a  function  of  depth  and  oxygen 
partial  pressure.  Additionally,  Figure  2.1  shows  acceptable 
mixture  ranges  for  a  given  depth  and  acceptable  depth  ranges 
for  a  given  mixture. 

Under  atmospheric  conditions,  the  partial  pressure  of 
oxygen  is  0.21  atm  since  the  percentage  of  oxygen  in  air  is 
approximately  21%.  Therefore,  at  the  surface,  divers  may 
safely  breathe  100%  oxygen  for  short  periods  (hours)  siuce  the 
oxygen  partial  pressure  is  only  1.0  atm,  well  belov/  the  1.6 
atm  upper  limit.  However,  as  depth  increases,  the  pressure 
also  increases  and  1.6  atm  can  be  quickly  exceeded  using  air 
as  the  diving-gas.  Equations  2.1  and  2.2  show  the 
relationship  between  pressure  (psia)  and  feet  of  sea  water 
(fsw)  . 
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OXYGEN  PERCENTAGE  vs.  DEPTH  (fsw) 


o  c 

•  '■  PERCENT  0^ 

( By  \'  o  1  u  n  1  e  ti  t  D  e  jT  i  l  i ) 
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DEPTH  (fsw) 

I'lLiurc'  2.\  Percentage  of  0^  in  diving-|[7as  mixture  as  a  function 

of  depth  and  oxygen  partial  pressure. 

(NOAA  Diving  Manual  1991) 
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P{psia) 


,  d{fsw)  ^-33 . 066  ^  (14,596) 
33.066 


2.1 


i 


» 


d(fsivr)  =  2.25P  -  33.066 


2.2 


Where, 

33.066  fsw  =  1  atm 
1  atm  =  14.696  psia 

Additionally,  the  partial  pressure  of  oxygen  for  given 
conditions  can  be  determined  from  equations  2.3  and  2.4. 


PPO2 


'  33 . 066 


+  1.0)  (%0.) 


2.3 


PPO2 


.  Pipsla) 
14.696 


)  (%02) 


2.4 


From  equation  2.3  it  can  be  shewn  that  the  greatest  depth  in 
which  air  can  be  safely  used  as  a  diving-gas  without  exceeding 
1.6  atm  partial  pressure  is  218  fsw. 

In  addition  to  quickly  exceeding  the  allowable  amount  of 
oxygen,  the  presence  of  nitrogen  renders  air  an  undesirable 
breathing  gas  at  depth.  Under  atmospheric  conditions  the 
partial  pressure  of  nitrogen  is  low  and  no  effects  are  felt. 
However,  when  the  partial  pressure  of  nitrogen  is  sufficiently 
high,  nitrogen  narcosis  or  rapture  of  the  deep  occurs. 
Nitrogen  narcosis  presents  symptoms  similar  to  alcohol 
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intoxication  and  can  be  extremely  dangerous  while  diving  at 
great  depths.  Table  2-1  indicates  the  events  which  may  occur 
while  diving  with  high  pressure  nitrogen  at  depth  (Randall 
1992)  .  Therefore,  air  is  not  a  suitable  diving-gas  and  is  not 
generally  used  for  depths  below  150  fsw.  The  advantage  of 
using  air  when  possible  is  that  it  is  inexpensive,  and  a 
chamber  can  be  ventilated  continuously  without  the  use  of 
carbon  dioxide  absorbers  which  are  addressed  ;  :r. 


Table 

2-1 

Effects 

of  Nitrogen  Narcosis  (Randall  1992) 

I# 

100 

-  150 

fsw 

light  head,  some  euphoria 

i 

150 

-  200 

fsw 

jovial,  some  dizziness 

1 

• 

200 

-  250 

fsw 

laughter  may  approach  hysteria 
simple  tasks  cannot  be  performed 
without  mistakes 

300 

fsw 

depression 

350 

fsw 

may  approach  unconsciousness 

• 

Since  helium  is  an  inert  gas  which  elindnates  the 
narcosis  effect,  it  is  often  used  at  depths  greater  than  150 
to  200  fsw.  However,  helium  does  have  the  side  effects  which 
are  listed  in  the  NOAA  Diving  Manual  (1991).  One  of  these 
side  effects  is  speech  distortion  which  increases  with  depth. 
Speech  distortion  can  be  overcome  by  using  helium  speech 
scramblers.  A  second  problem  with  helium  is  that  its  use 
results  in  rapid  body  heat  loss  since  the  thermal  conductivity 
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of  helium  is  six  times  that  of  air.  To  counter  this  problem, 
the  diving-gas  must  be  heated  for  dives  below  800  fsw.  A 
third  problem  associated  with  helium  is  the  high  pressure 
nervous  syndrome  (HPNS)  which  causes  tremors  in  the  hands  and 
jerky  movements  of  the  limbs.  This  problem  can  be  reduced  by 
using  a  slow  rate  of  compression.  Alternatives  to  helium  are 
also  listed  in  the  NOAA  Diving  Manual  (1991)  .  These 
alternatives  include  neon  and  hydrogen.  Neon  has  the 
advantages  over  helium  of  a  lower  thermal  conductivity  and 
less  speech  distortion.  However,  neon  is  much  more  expensive 
and  creates  more  breathing  resistance.  Hydrogen  is  easier  to 
breathe  at  depth,  but  it  causes  more  speech  distortion  than 
helium,  has  a  higher  thermal  conductivity,  and  can  be 
explosive.  Other  alternatives  may  be  under  development,  but 
are  beyond  the  scope  of  this  report. 

This  program  recommends  diving-gas  mixtures  based  upon  a 
helium-oxygen  mixture.  Smaller  percentages  of  oxygen  can  be 
mixed  with  helium  to  provide  a  diving-gas  mixture  which  can 
sustain  life  without  the  risk  of  oxygen  poisoning  at  deeper 
depths . 

Qaa  Mixtures 

The  ideal  diving-gas  mixture  is  determined  from  an  oxygen 
partial  pressure  of  0.35  atm.  Equation  2.5  calculates  the 
percent  oxygen  in  the  ideal  diving-gas  mixture. 


%0,  =  (  )  X  100% 

^  ■^ac/n 

Since  the  diving-gases  only  contain  helium  and  oxygen, 
^  equation  2.6  calculates  the  percent  helium  in  the  ideal 

diving-gas  mixture. 

•  =  100  -  %0, 


The  minimum  and  maximum  allowable  limits  for  oxygen  percentage 
are  calculated  by  equations  2.7  and  2.8  based  upon  oxygen 
partial  pressures  of  0.16  and  0.6  atm  respectively  as 
previously  discussed. 

n  16  2.7 

%o,  =  (iL-il)  X  100% 

^ac/n 


(.2-lI)  X  100% 

p 

^acm 


2.8 


If  the  maximum  allowable  oxygen  percentage  is  greater 
than  21%,  air  is  recommended  as  the  diving-gas  mixture. 
However,  in  most  cases  a  hyperbaric  chamber  operates  in  an 
environment  which  necessitates  the  use  of  a  helium-oxygen 
mixture. 

For  a  given  diving-gas  mixture,  the  pressure  (psia)  and 
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depth  (fsw)  ranges  are  calculated  using  equations  2.9  -  2.12. 


(14.696) 


(14.696) 


2.10 


<4in  =  (33.066)  -  33.066 


2.11 


t4ax  =  (33.066)  -  33.066 


2 . 12 


These  ranges  are  ai.i:;o  '-i  or/".  in  previous  figure  2.1. 


Ventilation  is  required  to  remove  carbon  dioxide  and 
replenish  oxygen  to  sustain  life.  When  air  is  used,  the 
required  ventilation  rate  must  be  sufficient  to  flush  the 
chamber  to  an  acceptable  level  of  carbon  dioxide.  This  rate 
shall  be  referred  to  as  the  "free  ventilation  rat.i".  The  free 
ventilation  rate  is  more  than  sufficient  to  provide  enough 
oxygen  to  sustain  life.  Since  air  is  relatively  inexpensive, 
the  free  ventilation  rate  is  desirable  when  utilizing  air  and 
carbon  dioxide  scrubbers  are  not  required  to  keep  the 
environment  carbon  dioxide  free. 
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The  free  ventilation  rate  must  be  sufficient  to  maintain 
an  environment  at  or  below  acceptable  partial  pressure  levels 
of  carbon  dioxide.  As  with  oxygen,  the  partial  pressure  of 
carbon  dioxide  increases  with  depth.  For  extended  periods  of 
time,  the  maximum  allowable  partial  pressure  of  carbon  dioxide 
should  not  exceed  0.005  atm.  Figure  2.2  shows  the  maximum 
allowable  CO2  partial  pressure  with  respect  to  mission 
duration.  It  is  evident  that  as  the  mission  duration 
increases,  the  smaller  the  carbon  dioxide  partial  pressure 
must  be  in  order  to  be  tolerated.  In  the  event  that  COj 
partial  pressure  exceeds  the  allowable,  hypercapnia  or  carbon 
dioxide  excess  can  jccur.  Table  2-2  lists  the  four  zones  of 
hypercapnia. 

Table  2-2  Effects  of  Hypercapnia 

Zone  I  -  No  effect  at  low  CO2  partial  pressures 

Zone  II  -  Minor  perceptive  changes 

Zone  III  -  Distracting  discomfort 

Zone  IV  -  Dizziness,  stupor,  unconsciousness 


Carbon  dioxide  production  is  determined  by  equation  2.13, 


CO2  =  (R^)  (0^)  (N^) 


2 . 13 
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EFFECTS  OF  CARBON  DIOXIDE 


► 


I 


CO.,  PARTIA:.,  PRESS!  RI-: 

(aim ) 


where , 


COjp  =  carbon  dioxide  production  (SCFM) 

Rq  =  respiratory  quotient  (0.9) 

Ojc  =  oxygen  consumption  (SCFM) 

Njj  =  number  of  divers 

Oxygen  consumption  is  determined  by  the  working  conditions 
listed  in  Table  2-3. 


Table  2-3  Oxygen  Consumption  of  Working  Conditions 


Condition 

Descriotion 

0..  consumotion 

1 

Resting 

0.5 

SLM 

/ 

0.0177 

SCFM 

2 

Light  Work 

0.7 

SLM 

/ 

0.0247 

SCFM 

3 

Moderate  Work 

1.5 

SLM 

/ 

0.0530 

SCFM 

4 

Heavy  Work 

2.0 

SLM 

/ 

0.0707 

SCFM 

5 

Severe  Work 

3.0 

SLM 

/ 

0.1060 

SCFM 

Figure  2.3  shows  the  oxygen  consumption  in  Standard  Liters  per 
Minute  (SLM)  for  each  of  the  working  conditions.  Standard 
conditions  exist  at  14.696  psia  and  70°  F. 

In  order  to  determine  the  free  ventilation  rate,  a  steedy 
state  condition  and  complete  mixing  are  assumed.  The 
volumetric  flow  rate  of  air  a^  depth  (V^j^^)  in  ft^/min  is 
determined  by  equation  2.14. 


( 


Ivc'lat  1  af  Respircitorv  Voluine  and  Oxygen  ('onsu. nipl  ion 
T\  pe  and  Levol  of  Rxeiiion  (Randall  1992) 


mRT 

PCO^ 


2.14 


where, 

m  =  the  mass  flow  rate  of  CO2  ( 

R  =  the  gas  constant  for  CO-  (35.11  Ib/lbmole-R) 

T  =  the  initial  temperature  (R) 

PCO2  =  the  partial  pressure  of  CO2  (psf) 

(0.005  atrn  or  10.58  psf  is  recommended  for  times 
greater  than  8  hours) 

The  ventilation  rate  (S^^^.)  in  Standard  Cubic  Feet  per  Minute 
(SCFM)  is  calculated  by  equation  2.15. 


14.696 


2.15 


The  mass  flow  rate  of  CO2  is  calculated  by  equation  2.16. 


(14.696psia)  (144psia/psf)  (co^) 

^  ^5.11  f  c~  lb/  lbmole-R)~  {53Q  R) 

where,  14.696  psia  and  530  R  are  standard  conditions 
Equation  2.16  can  be  simplified  to  equation  2.1'’. 

m  =  (0 . 1137)  (CO,  ) 


2. 16 


2  .  17 
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where  m  is  expressed  in  • 

Since  a  steady-state  condition  is  assumed,  it  is 
important  to  determine  now  long  it  takes  to  reach  this 

condition  at  the  beginning  of  a  mission  or  when  conditions 

change.  The  time  required  to  reach  a  95%  steady-state 

condition  can  be  derived  from  equation  2.18. 

PCO2  =  [1-e  ]  2.18 


where  -  the  interior  volume  of  the  chamber  (It^) 
Rearranging  to  solve  for  t  yields  equation  2.19. 


c 


V, 


ch 


V, 


In (1- ( 


V". 


PCO^ ) 


air 


mRT 


2.19 


By  substituting  in  the  values  for  PCO2  (0.005  atm) , 

R  (35.11  f t-lb/ Ibmole-R) ,  and  a  S5%  steady-state  condition, 
equation  2.19  becomes  equation  2.20. 


t  =  ( 


V. 


ch 


Un(: 


(.95)  (.005)  (14.7)  (144)  iV^ 
(35 . 11)  (/nT) 


)  ] 


2 .20 


Carrying  out  the  mult iplicatioii  yields  equation  2.21. 
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c  =  ( 


2 .21 


(-0.0167) 


V. 


mT 


where  the  factor  0.0167  converts  t  from  minutes  to  hours. 

I 

Carbon  Dioxide  Removal 

In  most  cases,  the  working  depth  or  pressure  requires  a 

I 

lielium-oxygen  mixture  for  reasons  previously  stated.  In  these 
cases  scrubbers  are  used  to  remove  carbon  dioxide  from  the 
chamber.  Therefore,  pure  oxygen  is  injected  but  care  must  be 
taken  to  ensure  proper  mixing  with  helium  in  the  chamber. 
This  mixing  occurs  while  the  gas  mixture  is  recirculating 
through  the  dehumidification  and  scrubbing  system. 

» 

The  amount  of  pure  oxygen  required  to  sustain  the  divers 
(02r)  if'  Standard  Cuoic  Feet  per  Minute  (SCFM)  is  determined 

from  equation  2.22. 

i 


Oy  =  ) 


(  *' 
V  n 


\ 

d' 


2 .22 


where  oxygen  consumption  (Oj^.)  was  previously  determined  from 
the  diver  working  conditions.  The  amount  of  CO2  produced 
(lbco2)  pounds  per  minute  is  calculated  by  equation  2.23. 


Ih^o,  " 


(60; 


2 .23 
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where  tjj  is  the  mission  duration  and  the  factor  of  60  converts 
minutes  to  hours. 

NAVFAC  DM-39  (1982)  lists  three  chemicals  used  as  carbon 
dioxide  absorbents.  These  chemicals  and  their  properties  are 
shown  on  Table  2-4.  Baralyme  is  no  longer  used  in  hyperbaric 
chambers  and  Lithium  Hydroxide  is  too  expensive  and  difficult 
to  handle  for  routine  use.  Therefore,  Sodasorb  is  the 
chemical  absorbent  used  in  this  program.  The  amount  of 
Sodasorb  (LBg)  required  for  each  mission  is  calculated  from 
equation  2.24  assuming  a  50%  efficiency. 


lb 


B 


0 . 245 


2.24 


Although  it  is  possible  to  obtain  higher  absorbent 
efficiencies,  canister  design,  temperature,  humidity,  and 
moisture  condensation  must  be  considered.  The  50%  efficiency 
is  a  conservative  estimate  used  for  design  purposes.  During 
actual  chamber  operations,  chemical  detection  tubes  must  be 
used  to  accurately  determine  when  to  change  carbon  dioxide 
absorbent  canisters.  The  number  of  carbon  dioxide  absorbent 
canisters  required  for  a  mission  will  depend  on  the 
capabilities  of  individual  hyperbaric  chambers. 


Dlvina-Qas  Mixture  Properties 

Properties  of  the  diving-gas  mixture  must  be  accurately 
determined  in  order  to  calculate  energy  and  dehumidification 
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Table  2-4  Characteristics  of  Carbon  Dioxide  Absorbents 


(NAVFAC  DM- 3 9 

1982) 

Absorbent 

Characteristic 

Baralyme 

Lithium 

hvdroxide 

Sodasorb 

Density  (Ib/ft^) 

65.4 

28.0 

55.4 

Theoretical  COj 
absorption  (lb  COj/lb) 

0.39 

0.92 

0.49 

Useful  CO2  absorption  0.195 

(based  on  50%  efficiency) 

(lb  COj/lb) 

0.46 

0.245 

requirements  which  will  be  discussed  later.  Properties  of 
concern  include  molecular  weight,  specific  heat,  and  density 
of  the  selected  diving-gas  mixture. 

The  molecular  weight  of  a  substance  is  defined  as  the 
weight  in  pounds  of  one  pound-mole  of  the  substance.  The 
molecular  weight  of  the  diving-gas  mixtures  is  dependent  upon 
the  percentages  of  oxygen  and  helium  in  the  mixture  and 
therefore,  decreases  as  the  percent  oxygen  in  the  mixture 
decreases.  This  in  turn  decreases  the  density  which  minimizes 
the  energy  requirements.  If  the  diving-gas  is  air,  the 
molecular  weight  is  29.0  Ib/lb-mole.  Otherwise,  the  mole 
fractions  of  oxygen  and  helium  in  the  diving-gas  mixture  must 
first  be  determined.  Mole  fractions  (X)  are  determined  by 
equations  2.25  and  2.26.  The  sum  of  the  mole  fractions  must 
always  equal  1, 


2.26 


Once  mole  fractions  are  determined,  the  molecular  weight  in 
lb/ lb-mole  can  be  calculated  from  equation  2.27. 


m  =  [  {Xq)  {mw^)  ]  +  [  (x„^)  (Mw„j  ] 


where, 

MWo2  =  32.00  lb/ lb-mole 
MWHe  =  4.0026  Ib/lb-mole 

The  specific  heat  at  constant  pressure  (Cp) ,  in  units  of 
BTU/ lb-mole,  of  the  diving-gas  mixture  must  be  known  for 
determining  energy  requirements.  Cp  is  defined  by  equation 
2.28. 


C  =  (  —  ) 
'-p 


2.28 


where,  h  is  the  enthalpy  and  T  is  temperature. 

Given  the  mole  fractions  previously  determined,  Cp  in  BTU/ lb- 
mole  can  be  calculated  from  equation  2.29. 


2 . 29 
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where , 

Cpo2  =  6.95  BTU/ lb-mole 
CpHe  =  4.965  BTU/ lb-mole 

When  air  is  used  as  the  diving-gas,  Cp  is  6.95  BTU/ lb-mole. 
Table  2-5  lists  the  molecular  weights  and  specific  heats  of 
common  helium/oxygen  mixtures  and  air. 


a 

Table  2-5  Properties 

of  Common  Diving-Gas 

Mixtures 

w 

Diving-Gas  Mixture 
(%OW%He) 

(BTU/lb^ole) 

MW 

f lb/ lb-mole^ 

1/99 

4-985 

4.283 

5/95 

5.064 

5.403 

10/90 

5.164 

6.802 

15/85 

5.263 

8.202 

20/80 

5.362 

9.602 

40/60 

5.759 

15.202 

Air 

6.950 

29.000 

• 

Density  of  air  is  calculated  using  the  ideal  gas  equation 
of  state  since  air  will  only  be  used  down  to  150  fsw  (81.36 
psia) .  The  ideal  equation  of  state  is  shown  as  equation  2.3L. 


PV  =  nR^T 


2 . 30 


25 


where, 

P  =  absolute  pressure  (psf) 

V  =  volume  (cu  ft) 
n  =  number  of  moles 

Ry  =  Universal  Gas  Constant  (1544.0  ft-lb/lbmole-R) 

T  =  absolute  temperature  (R) 

Equation  2.30  can  be  simplified  and  rearranged  to  solve  for 
density  in  equation  2.31. 

./Trt  /  2.31 

_  kMW)  (P) 

R,T 

Density  of  helium-oxygen  mixtures  is  determined  from  the 
virial  equation  of  state.  Dymond  and  Smith  (1980)  present  the 
virial  equation  of  state  shown  as  equation  2.32. 


R,T 


=  1  + 


S(D 


C(T) 


D{T) 


2.32 


where,  the  values  of  P/Ry,  and  T  are  the  same  as  those  used 
above,  and 

=  molar  volume  (cu  ft/lb-mole) 

1  =  the  first  virial  coefficient 
B(T)  =  the  second  virial  coefficient 
C(T)  =  the  third  virial  coefficient 
D(T)  =  the  fourth  virial  coefficient 


# 


It  is  common  practice  to  calculate  the  density  of  gases  based 
upon  only  the  first  three  virial  coefficients  since  data 
beyond  this  range  are  scarce.  Dymond  and  Smith  (1980)  state 
that  uncertainties  in  the  second  virial  coefficient  are 
relatively  small  (usually  less  than  1%)  but  are  much  larger 
for  the  third  virial  coefficient.  Additionally,  the  third 
virial  coefficient  is  much  more  temperature  dependent. 

Dymond  and  Smith  (1980)  developed  the  most  recent  table 
for  determining  the  second  and  third  virial  coefficients  for 
gas  mixtures.  Equation  2.33  calculates  the  second  virial 
coefficient  of  a  gas  mixture. 


5„(T)  =  B. 

fn  '  X 


(T)  (X.)^  f  B.,(r)  (X.) 


x^) 


f  B„  (D  (X,)- 


2.33 


where, 

B^(T)  =  the  second  virial  coefficient  of  the  mixture 
Xj^  =  the  mole  fraction  of  the  first  pure  gas 
X2  =  the  mole  fraction  of  the  second  pure  gas 

=  the  second  virial  coefficient  of  the  first  pure  gas 
B22  =  the  second  virial  coefficient  of  the  second  pure  gas 
®12  ~  mixed  virial  coefficient 
Equation  2.34  calculates  the  third  virial  coefficient  of  a  gas 
mixture. 


27 


i 


=  ^111  (2’)  ^  3C,,2(r)  (Xi)2(xj 

^  3C,22(T)  (A-,)  (X^)^  -H  C^^^iT)  ix.y 


2.34 


# 


# 


where, 

C^(T)  =  the  third  virial  coefficient  of  the  mixture 

-  the  third  virial  coefficient  of  the  first  pure  gas 
Cii2  =  the  contribution  of  two  molecules  of  the  first  pure 
gas  and  one  molecule  of  the  second  pure  gas 
^122  ~  contribution  of  one  molecule  of  the  first  pure 
gas  and  two  molecules  of  the  second  pure  gas 
Due  to  the  absence  of  data  for  determining  the  third  virial 
coefficients  for  helium-oxygen  mixtures,  the  density  is 
calculated  based  upon  the  first  two  virial  coefficients  only. 
Tables  2-6  lists  the  virial  coefficients  for  oxygen,  helium, 
and  a  helium-oxygen  mixture  respectively  from  Dymond  and  Smith 
(1980).  Since  this  is  the  most  recent  data  available,  it  is 
obvious  that  very  little  data  exists  to  determine  the  mixed 
virial  coefficient  for  helium-oxygen  mixtures. 

By  substituting  the  values  from  Tables  2-6  into  equation 
2.33,  the  second  virial  coefficient  in  SI  units  (cm^/g-mole) 
is  calculated  from  equation  2.35. 


5=  [  (11 .44)  -  [  (8.8)  (A'„  )  ) 


[  (16.24)  (Ay^)"] 


2 . 35 
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Table  2-6  Virial  Coefficients  of  Helium,  Oxygen,  and 
Heliox  Mixtures 
(Dymond  and  Smith  1980) 


Virial  Coefficients 

Helium  Oxygen  Mixture 


TfC) 

B 

C 

B 

C 

B 

C 

-183 

10.76 

277 

-241 

-12764 

-4.4 

- 

0 

11.87 

75.5 

-21.89 

1230 

- 

- 

25 

11.44 

72.1 

-16.24 

1163 

- 

- 

50 

11.37 

72.3 

-11.62 

1146 

- 

- 

B  = 

cm^/g-mole  C 

=  cm^/g-mole^ 

In  order  to  use  this  coefficient,  it  must  be  converted  into 
the  English  units  of  ft^/ lb-mole.  This  conversion  is  made  by 
equation  2.36. 


B  =  ( 


ifc- 


)(^) 


28316. 85c/n’  lb 

=  (0.0160)  {Bsj) 


2 .36 


Solving  the  virial  equation  of  state  for  molar  volume  (ft^/ lb- 
mole)  yields  equation  2.37. 


(K  «  MK£? 
2 


2  .  3  7 
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where  K  is  a  constant  defined  as  R^T/P. 

Finally,  density  is  calculated  by  equation  2.38. 


P 


m 


2 . 38 


where  MW  was  calculated  in  equation  2.27  and  molar  volume  was 
calculated  in  equation  2.37. 

The  U.  S.  Navy  Diving  Gas  Manual  (1971)  lists  the 
densities  of  standard  diving-aa*^  mixtures  at  varying 
conditions.  The  resulting  err'jr  irom  these  equations  is  less 
than  0.2%. 

Psvchromatric  Charts 

Psychrometric  charts  are  used  to  compute  heating  and 
cooling  requirements  along  with  dehuroidif ication  and  moisture 
removal  requirements.  Standard  psychrometric  charts  can  be 
used  for  calculating  the  effects  of  heating  and  cooling  moist 
air  at  atmospheric  conditions  but  cannot  be  used  with  diving- 
gas  mixtures  at  much  greater  depths.  Although  these  charts 
are  aval  ...able  for  pressures  up  to  600  psia  (approx.  1300  fsw) 
in  the  U.  S.  Navy  Diving  Gas  Manual  (1971),  they  are  time 
consuming  and  lack  the  accuracy  of  this  computer  program. 
Therefore,  it  is  necessary  to  develop  relationships  which  can 
be  used  in  computer  simulation.  Additionally,  reading  values 
off  these  charts  requires  non-linear  interpolation  between  100 
psia  curves.  Figure  2.4  shows  the  psychrometric  chart  in  the 
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U.  S.  Navy  Diving  Gas  Manual  (1971)  developed  for  gas  mixtures 
at  pressures  from  100  to  600  psia. 

Given  the  dry  bulb  temperature  and  saturation  lines  at 
constant  pressure,  the  moisture  mole  fraction  (moles  water 
vapor /moles  dry  gas)  can  be  determined.  The  dry  bulb 
temperature  is  defined  by  ASHRAE  (1967)  as  the  Fahrenheit 
temperature  of  moist  gas  at  rest  with  respect  to  the 
temperature-measuring  element.  The  saturation  lines  at 
constant  pressure  represent  100%  relative  humidity  at  that 
pre.ssure . 

The  moisture  mole  fraction  in  moles  water  vapor/mole  dry 
gas  is  determined  by  equation  2.39  if  the  pressure  of  the 
environment  and  saturation  pressure  are  known. 


2.39 


where, 

=  the  mole  fraction  for  water  vapor  at  saturation 
(moles  water  vapor/mole  dry  gas) 

Pg  =  the  vapor  pressure  of  water  in  saturated  moist  air 
(psia) 

P  =  the  absolute  pressure  of  the  chamber  (psia) 

The  values  of  Pg  have  been  tabulated  by  ASHRAE  (1967)  for 
various  temperatures.  These  values  are  reproduced  in  Table 
2-7.  Equation  2.40  is  a  best  fit  curve  through  these  points 
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Table  3-7  Vapor  Pressure  of  Water  in  Saturated  Moi>_t  Air  (psia) 
ASHRAE  (1967) 


I 


i 


Absolute 

Pressure  (Ps) 

1 

Temp  (F) 

;  Ps  (psia) 

1 

1 

i  Temp  (r) 

1 

1 

Ps  (psia)  ■ 

1 

40 

1 

1 

1  0.12164 

1 

!  72 

I  ; 

1  ! 

1  ; 

0.38356  ' 

41 

j  0.12646 

73 

0.40190 

42 

0.13145 

i 

0  41564  I 

43 

1  0.13660 

'■  75 

1  0.42979  j 

44 

:  0.14  194 

76 

■  0  44  435  1 

45 

i  o:  14746 

■  77 

0.45935  ; 

46 

1  0  15317 

78 

i  0.47478  : 

47 

1  0  15907 

:  79 

0.49066 

48 

016517 

:  80 

0  50701 

;  -^9 

0.17148 

81 

0.52382 

i  50 

1  0.17799 

j  82 

j  0.54112  ; 

;  51 

1  t  q 

1  w  i  J  •  r  o 

83 

0.55892  ; 

i  52 

0  19169 

84 

1  0.57722  : 

53 

;  0.19888 

1  85 

1  0.59604 

54 

j  0.20630 

'  86 

'  0.61540  I 

55 

!  0.21397 

87 

0  63530  ■ 

56 

0  22188 

;  88 

i  0  65575 

i  57 

0.23006 

;  89 

0.67678 

1  58 

0.23849 

90 

:  0.69838 

59 

0  24720 

;  91 

0,72059 

■  60 

r\  I'\  —  r'  ^  !■» 

U-40U  1  () 

:  92 

0  74340 

,  61 

:  0.26545 

93 

1  0  76684 

62 

0,27502 

94 

!  0.79091 

63 

0  28488 

9:5 

0  81564 

■  64 

0  29505 

9(i 

0  84  103 

65 

0  30554 

97 

0  86711 

i  66 

1  031636 

98 

0  893418 

67 

0.:'12750 

99 

0  92  617 

68 

0  33900 

100 

()  9.i959 

69 

0  35084 

1  70 

036304 

71 

■  037561 

determined  from  PCMATLAB  (1989)  with  a  correlation  coefficient 
of  0.9999S65. 

Pg  =  (1.9285  X  lO'^M  T®  +  (9.8685  x  10-^°)  T* 

+  (3.3193  X  10-7)  +  (1.9916  x  lO'^)  2.40 

+  (1.04  88  X  10'^)  r-^  (2.2080  xlO’2) 

Where  T  =  the  chamber  temperature  {°  7)  . 

The  U.  S.  Navy  Diving  Gas  Manual  (1971)  states  that  for 
total  mixture  pressures  up  to  100  psia,  no  correctio..  for  the 
effect  of  total  pressure  on  the  water-vapor  saturacion 
pressure  is  used.  However,  above  100  psia,  Poynting's 
equation  is  used  as  an  estimated  correction  for  total 
pressure.  Casstone  (1947)  derives  Poynting's  Equation  shown 
as  equation  2.41. 

_  V-ilAP)  (P.) 

^  p„r  2.41 


where  is  the  specific  volume  of  the  saturated  liquid. 
Values  of  Vj^  are  tabulated  by  ASHRAE  (1967)  and  are  reproduced 
in  Table  2-8  in  units  of  (ft^/lb) .  To  be  used  in  Poynting's 
Equation,  must  be  expressed  in  (ft^/mole)  .  This  is  easily 
accomplished  by  multiplying  Vj^  by  the  molecular  weight  of 
water  (18.0158  Ib/lb-mole).  Then,  the  corrected  saturation 
pressure  will  be  calculated  from  equation  2.42. 
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Table  2-8  Specific  Volume  of  Saturated  Liquid 
ASHRAE  (1967) 


Specific  Volume  (cu  ft,  lb) 


Temp  (F)  ! 

1 

1 

1 

V  i 

Temp  (F) 

V, 

40 

0  01602  ! 

72 

0  01606 

41  I 

0.01602  j 

73 

0.01606 

42 

0.01602  ! 

74 

0  01606 

43  1 

0.01602  1 

75 

0  01606 

44  ; 

0.01602  ; 

76 

0.01606 

45  ! 

0.01602  1 

77 

0.01607 

46  : 

0.01602  i 

78 

0.01607 

47 

0  01602  ' 

79 

0.01607 

48  ; 

0,01602 

80 

0  0  1607 

49  ; 

0.01602  ; 

81 

0  01608 

50  : 

0.01602  ! 

82 

0  01608 

51  ' 

0.01602 

83 

0.01608 

52  i 

0.01602  ; 

84 

0,01608 

53  ! 

0.01603  ; 

85 

0,01609 

54  1 

0.01603  i 

86 

0.01609 

55  1 

0  01603  ■ 

87 

C. 01609 

56  ; 

0  01603  : 

88 

0.016'^ 

57  ! 

0.01603 

89 

0.016 

53 

0  01603 

90 

0  0  16K. 

59 

O  n  1  > 

n  1 

0  01610 

60 

0  01603 

92 

0.0  161  ! 

6)  ; 

0.01604  1 

93 

0.0161  1 

62  ; 

001604 

9.) 

0  0161  1 

63 

0  01604 

95 

0  0  1  6  1 

64 

0  01604  : 

96 

0  0  1612 

65 

0.01604 

97 

001612 

66 

0  01601 

98 

0  0  161.2 

67 

0  016) 5 

99 

0  0  6  i  :  1 

68 

0  0 1605 

100 

0  01613 

69 

0  0 1605 

70 

0  0 1 605 

71 

001605 

2.42 


► 


Pa  =  Pa 


AP. 


Where  Pg^  is  the  uncorrected  saturation  pressure.  The 
moisture  mole  fraction  of  water  vapor  at  saturation  (X^g)  is 
then  calculated  from  equation  2.39. 

The  molar  enthalpy  of  moisture  content  at  saturation  (H) 
in  (BTU/mole  dry  gas)  is  calculated  from  equation  2.43. 


» 


H  =■-  (20, 000) 


2.43 


I 


Enthalpy  is  a  thermodynamic  property  defined  by  equation  2.44. 


H  =  U  ^  PV 


2.44 


where, 

U  =  internal  energy 
P  =  absolute  pressure 
V  =  volume 

Determination  of  enthalpy  is  important  for  calculating  energy 
requirements. 

Dehumidif ioation  and  Energy  Requirements 

In  order  to  determine  dehui,.xdif ication  requirements,  the 
initial  and  desired  relative  humidities  must  be  known.  Given 
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the  initial  and  desired  relative  humidities,  the  moisture  mole 
fraction  and  molar  enthalpy  of  moisture  content  can  be 
determined  for  both  conditions  by  equations  2.45  -  2.48. 


where  RHj^  and  are  the  initial  and  desired  relative 

humidities  expressed  in  fractions,  and  X^g  and  H  were 
previously  determined.  Once  the  moisture  mole  fraction  at  the 
desired  relative  humidity  is  known,  the  maximum  coil 
temperature  required  for  dehumidification  can  be  calculated. 
Equation  2.49  calculates  the  new  saturation  pressure  (Pg) . 

(X^  )  (P) 

P  =  — l:^ - 

®  1  +  2.49 
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Then,  the  maximum  coil  temperature  (T^.)  in  (°F)  can  be 
determined  as  the  intersection  of  saturation  pressure  curve  at 
X^2  from  equation  2.50. 


=  (2 .4135  X  10^)  P/  -  (8.4538  X  102)  P/ 

+  (1.1635  X  10^)  P,^  -  (8.1910  X10^)P/  2.50 

+  (3 . 5407  X  10^)  P,  (7.0106) 


Equation  2.50  was  calculated  using  PCMATLAB  (1989)  and  has  a 
correlation  coefficient  of  0.999974. 

Once  the  coil  temperature  is  determined,  the  moist  gas 
mixture  enthalpies  at  the  initial  and  coil  temperatures  are 
calculated.  Equations  2.51  and  2.52  calculate  the  moist  gas 
mixture  enthalpies  in  (BTU/lb)  at  the  initial  and  coil 
temperatures  respectively. 


[(cp(rj  *//,] 

MW 


2.52 


where  the  specific  heat  and  molecular  weight  were  previously 
determined  as  properties  of  the  diving-gas  mixture. 

The  energy  (Q^)  in  (BTU/min)  required  to  dehumidify  the 
chamber  to  the  desired  relative  humidity  is  calculated  from 
equation  2.53. 
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where,  is  either  the  free  ventilation  rate  for  air 
systems  or  the  recirculation  flowrate  for  helium-oxygen 
systems.  Since  the  coil  temperature  is  generally  colder  than 
the  desired  diving-gas  mixture,  the  diving-gas  must  be 
reheated  to  a  desired  temperature.  The  change  in  enthalpy 
from  the  coil  to  a  desired  temperature  (T^j)  is  calculated  by 
equation  2.54. 


Ah 


l(C^)  -  T^)] 

MW 


2.54 


Then,  the  energy  required  to  reheat  the  diving-gas  (Q^.)  to  the 
desired  temperature  is  calculated  by  equation  2.55. 


0.  = 


2.55 


The  total  energy  (Q^)  required  by  the  system  for 
dehumidification  and  reheat  is  calculated  by  equation  2.56. 


2 . 56 


However,  energy  lost  to  heat  transfer  through  the  chamber 


walls  depends  on  insulation  and  is  not  considered  in  this 
analysis.  It  would  be  an  additional  energy  requirement. 

Moisture  Removal  Rate 

In  all  dehumidif ication  processes,  moisture  removal  must 
be  addressed.  Due  to  the  rapid  cooling  of  the  moist  air, 
water  vapor  condenses  and  must  be  removed  from  the  chamber. 
In  order  to  determine  the  moisture  removal  rate,  the  humidity 
ratios  at  the  initial  and  desired  relative  humidities  must  be 
determined.  Equations  2.57  and  2.58  calculate  the  humidity 
ratios  in  pounds  water  vapor  per  pounds  dry  air  (Ib^^/lb^j)  at 
the  initial  (W^)  and  desired  (Wj)  relative  humidities 
respectively. 


(A'  )  (18.0158) 
MW 


2.57 


{X^)  (18.0158) 
MW 


2.58 


where  18.0158  is  the  molecular  weight  of  water.  The  moisture 
removal  rate  (MRR)  in  (Ib/day)  can  then  calculated  by  equation 
2.59. 


MRR  =  (A^^  -  (Vj)  V'^^^p(i440)  lb! day 


2.59 


where, 


is  either  the  free  ventilation  rate  or  the 


recirculation  flowrate  and  the  factor  of  1440  converts 
Ib/minute  to  Ib/day.  This  moisture  removal  rate  results  from 
the  dehuroidif ication  process  only.  The  moisture  produced  by 
the  carbon  dioxide  absorbents  is  not  included  in  this 
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DESCRIPTION  OP  SIMULATION 

Simulation  Structure 

This  program  is  designed  to  simulate  hyperbaric  chamber 
environments  such  as  decompression  chambers,  submarine  lockout 
and  pilot  chambers,  diving  bells,  underwater  habitats,  and 
ocean  simulation  facilities  under  different  conditions.  The 
theory  and  equations  developed  in  Chapter  II  are  incorporated 
into  this  interactive  computer  program  which  is  written  in 
Fortran  77  for  a  386  personal  computer.  In  order  to  design  a 
hyperbaric  chamber  with  this  program,  the  input  listed  in 
Table  3-1  is  required.  Given  this  input,  the  program 
calculates  values  and  provides  the  output  listed  in  Table  3-2. 

A  flowchart  for  this  program  is  included  as  Figure  3 . l 
and  the  source  code  is  provided  as  Appendix  B. 

gubroutlne  Desoriptlop 

This  program  makes  use  of  several  subroutines  to 
calculate  the  required  quantities.  A  description  of  each 
subroutine  follows  in  the  paragraphs  below. 

Subroutine  ’’CONVERT"  converts  depth  to  pressure.  For  all 
subsequent  calculations,  this  program  uses  absolute  pressure. 
If  the  user  chooses  to  input  a  working  depth  in  feet  of 
seawater  (fsw),  "CONVERT"  converts  depth  to  pressure  using 
equation  2.1.  This  value  is  returned  to  the  main  program  and 
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Coll  Tenipcr«iurr 


Figure  3  ]  Flowchart  for  "Simulation  of  a  Hyperbaric- 
Chamber  Environment" 


Table  3-1  Required  Program  Input 


Reouirement 

Units 

1) 

Working  Depth  fsw  (feet 

or  Pressure 

of  seawater) 
psia 

2) 

Initial  Temperature 

0  p 

3) 

Desired  Temperature 

0  p 

4) 

Initial  Relative  Humidity 

% 

5) 

Desired  Relative  Humidity 

% 

6} 

Number  of  Divers 

N/A 

7) 

Internal  Volume  of  Chamber 

ft3 

8) 

Diver  Working  Conditions 

N/A 

9)* 

Mission  Duration 

hrs 

10) 

Recirculation  Flowrate 

ft^/min 

*  -  These  values  are  only  required  if  CO2 
going  to  be  used  during  the  mission. 

scrubbers  are 

used  for  all  future  calculations. 

Subroutine  "OASMIX”  calculates  the  minimum  and  maximum  oxygen 
percentages  for  the  pressure  of  interest  based  upon  oxygen 
partial  pressures  of  0.16  and  0.6  atm  respectively  from 
equations  2.7  and  2.8.  If  the  maximum  acceptable  oxygen 
percentage  is  greater  than  21%,  the  ideal  diving-gas  mixture 
is  air.  All  further  calculations  are  then  based  upon  air  as 
the  diving-gas  mixture.  If  the  maximum  acceptable  oxygen 
percentage  is  less  than  21%,  the  program  determines  an 
acceptable  oxygen  percentage  using  an  IF-THEN  statement. 


Table  3-2  Program  Output 


OUtDUt 

Pnits 

1) 

Diving-Gas  Mixture 

%He/%02 

2) 

Oxygen  Required 
or  Ventilation  Rate 

SCFM 

3*) 

Sodasorb  Required 

lbs 

4) 

Time  to  Reach  95%  Steady-State 

hrs 

5) 

Molecular  Weight 

lb/ lb-mole 

6) 

Gas  Mixture  Density 

Ib/ft^ 

7) 

Maximum  Coil  Temperature 

O  p 

8) 

Dehumidification  Energy  Requirements 

BTU/min 

9) 

Reheat  Energy  Requirements 

BTU/min 

10) 

Moisture  Removal  Rate 

Ib/day 

-  This  value  is  only  provided  if  COj 
used  during  the  mission. 

scrubbers  are 

Siibroutine  "OASRANOE"  calculates  a  pressure  and  depth  range 
based  upon  equations  2.9  -  2.12  after  an  acceptable  diving-gas 
mixture  has  been  selected.  This  information  is  useful  if  the 
same  diving-gas  mixture  is  used  for  multiple  depths.  If  the 
user  is  planning  to  operate  outside  these  limits,  a  different 
diving-gas  mixture  should  be  selected. 

Subroutine  "PREERATE"  determines  the  ventilation  rate  required 
to  keep  the  carbon  dioxide  partial  pressure  at  an  acceptable 
level  when  air  is  used  as  the  diving-gas.  First,  "FREERATE" 
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calculates  the  carbon  dioxide  production  based  upon  the 
working  conditions  and  equation  2.13.  In  order  to  maintain  a 
CO2  partial  pressure  below  0.005  atm,  the  volumetric  flow 
rates  of  air  at  depth  and  at  the  surface  are  calculated  from 
equations  2.14  and  2.15. 

Siibroutine  ''AB8RATE'*  calculates  the  amount  of  oxygen  required 
to  sustain  life,  the  amount  of  CO2  produced,  and  the  amount  of 
Sodasorb  required  to  maintain  an  acceptable  environment  when 
the  diving-gas  is  a  helium-oxygen  mixture.  The  amount  of 
oxygen  required  is  calculated  by  equation  2.22  while  the 
amount  of  CO2  produced  in  pounds  is  calculated  by  equation 
2.23  and  the  amount  of  Sodasorb  required  is  calculated  by 
equation  2.24.  The  user  may  wish  to  design  different  CO2 
scrubbing  systems  based  upon  the  amount  of  CO2  produced. 

Subroutine  "STSTATE'*  calculates  the  time  required  to  reach  a 
95%  steady-state  condition.  Since  steady-state  cond ;  *-.'.or5  are 
assumed,  it  is  important  to  determine  the  time  reqxiired  to 
reach  this  condition.  The  time  required  to  n^ach  the  95% 
steady-state  condition  is  calculated  by  equation  2.21.  This 
information  is  useful  when  beginning  a  mission  or  changing 
conditions . 

Subroutines  '*M0LWT",  "SPECHEAT”,  and  ’’DENSITY"  calculate  che 
properties  of  the  selected  diving-gas  mixture.  If  the  diving- 
gas  mixture  is  air,  the  molecular  weight  and  specific  heat  are 


constant  values  of  29.0  Ib/lb-mole  and  6.95  BTU/ lb-mole 
respectively. 

However,  for  helium-oxygen  mixtures,  "MOLWT"  calculates 
the  molecular  weight  from  equation  2.27  and  "SPECHEAT*' 
calculates  the  specific  heat  from  equation  2.29.  This 
information  is  used  in  subsequent  calculations. 

'•DENSITY"  calculates  the  density  of  the  diving-gas  for 
both  air  and  helium-oxygen  mixtures.  Since  air  is  only  used 
at  low  pressures,  "DENSITY"  determines  the  density  of  air  from 
the  ideal  equation  of  state,  equation  2.31.  When  helium- 
oxygen  mixtures  are  used,  "DENSITY"  calculates  density  from 
the  virial  equation  of  state,  equation  2.32.  The  second 
virial  coefficient  for  the  mixture  is  determined  from  equation 
2.35  and  converted  to  English  units  in  equation  2.36.  This 
value  is  then  used  to  determine  the  molar  volume  from  equation 
2.37,  and  finally  the  density  is  calculated  by  equation  2.38. 
This  value  is  used  in  subsequent  energy  calculations. 

Subroutine  "PSYCHART"  calculates  the  moisture  mole  fraction 
(X^^g)  and  the  enthalpy  of  moisture  content  (H)  for  the  diving- 
gas  mixture.  First,  the  vapor  pressure  of  water  in  saturated 
moist  air  (Pg)  is  calculated  from  equation  2.40.  If  the 
working  pressure  is  less  than  100  psia,  no  correction  is 
required.  However,  for  pressures  above  100  psia,  Pg  is 
corrected  using  Poynting's  equation  2.41  where  the  specific 
volume  of  the  saturated  liquid  (V^^)  is  read  off  of  Table  2-8 
using  an  IF-THEN  statement.  The  corrected  value  of  Pg  is  then 


calculated  by  equation  2.12.  The  values  of  and  H  are 
calculated  from  equations  2.39  and  2.43  respectively. 

Subroutine  ’’XWSH"  calcv.iaJ'es  the  initial  and  desired  moisture 
mole  fractions  (X^j^  and  X^.,)  and  molar  enthalpies  and  Hj) . 
These  properties  are  determined  by  the  initial  and  desired 
relative  humidities.  The  values  are  calculated  based  upon  X^g 
and  H  by  equations  2.45  -  2.48. 

Subroutine  "COILTEMP”  calculates  the  maximum  coil  temperature 
required  to  dehumidify  the  chamber  to  the  desired  relative 
humidity.  The  moisture  mole  fraction  at  the  desired  relati/e 
humidity  (X^^2)  calculated  by  equation  2.46.  The  new 
saturation  pressure  (Pg)  is  determined  by  solving  equation 
2.39  for  Pg  which  leads  to  equation  2.49.  Given  Pg,  the  coil 
temperature  (T^,)  is  calculated  by  equation  2.50. 

subroutine  "MaMENTH"  calculates  the  moist  gas  mixture 
enthalpies  entering  the  coolirg  coil  (ll^)  and  leaving  the 
cooling  coil  (H^)  by  equationu  2.51  and  2.52  respectively.  In 
this  program,  the  temperature  entering  the  cooling  coil  is 
as  umed  to  be  the  initial  temperature  (Tj^)  and  the  temperature 
leaving  the  cooling  coil  is  assumed  to  be  the  maximum  coil 
temperature  (T^,)  .  These  values  may  differ  in  an  actua] 
cooling  system. 

Subroutine  "ENERGY"  calculates  the  energy  required  lor 
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dehumidif ication  and  reheat  *:<.  desired  temperature.  The 
energy  required  for  dehumidif ' ^at ion  (Q^j)  is  calculated  by 
equation  2.53.  If  the  diving-gas  is  air,  the  free  ventilation 
rate  calculated  by  "FREERATE”  is  used.  For  helium- 
oxygen  mixtures,  the  recirculation  flowrate  input  by  the  user 
is  used  as  In  order  to  determine  the  energy  required  to 
reheat  to  a  desired  temperature  (Qj.) ,  the  change  in  enthalpy 
from  the  coil  temperature  to  a  desired  temperature  must  be 
calculated  by  equation  2.54.  Then,  Q^.  is  calculated  in  a 
similar  manner  as  by  equation  2.55.  The  total  energy 
required  (Q^)  is  calculated  by  equation  2.56. 

8\ibroutin«  "MOISTREM"  caici.lates  the  moisture  removal  rate 
required  for  the  above  dehumidif ication.  First,  the  humidity 
ratios  at  the  initial  (W^)  and  desired  (W2)  relative 
humidities  are  calculated  by  equations  2.57  and  2.58.  Then, 
the  moisture  removal  rate  is  calculated  by  equation  2.59  where 
again,  is  the  free  ventilation  rate  for  air  or  the 
recirculation  flowrate  for  helium-oxygen  mixtures.  The 
calculations  performed  by  these  subroutines  are  incorporated 
into  the  main  program,  and  the  results  are  displayed  to  the 
user  as  program  "output”.  Three  examples  are  discussed  in 
detail  in  the  next  chapter. 


MiHlIMlHIIIIHI 
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CHAPTER  IV 


APPLICATION  OF  SIMULATION 

This  program  can  be  used  not  only  to  simulate  a 
hyperbaric  chamber  but  also  a  submarine  diver  lock-out  chamber 
or  an  underwater  habitat.  The  following  examples  demonstrate 
this  program's  usefulness  for  determining  the  requirements  of 
these  cases. 


Exanple  l: _ Hyperbaric  chanber 

This  example  simulates  a  hyperbaric  chamber  used  for  a 
deep  dive.  As  a  model,  the  tJ.  S.  Navy's  Ocean  Simulation 
Facility  at  the  Naval  Coastal  Systems  Center  in  Panama  City, 
FI  was  be  used.  This  facility  as  shown  in  Figure  4.1  is 
capable  of  pressurization  down  to  2250  fsw  and  is  described  in 
detail  by  Shilling  et  al.  (1976)  and  U.  S.  Navy  Diving  Manual 
(1991) .  The  input  conditions  for  this  example  are  given  in 
Table  4-1.  The  calculation  steps  are  demonstrated  below: 

1)  Since  a  working  depth  was  input,  it  was  converted  to 
a  pressure.  By  equation  2.1, 


P  = 


.  (2000  ^  33.066) 
33.066 


) (14.696) 


P  =  903.6  psia 
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Figure  4.1  Ocean  Simulation  Facility 
(Shilling  et  al.  1976) 


2)  The  ideal  gas  mixture  for  this  depth  was  calculated 
by  equations  2.5  and  2.6. 


%02 

%02 


( 


0.35 


( 


903.6 


14 
0.57% 


.696 


) (100%) 


%He  *  100  -  0.57 
%He  =  99.43% 
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Table  4-1  Hyperbaric  Chamber  Simulation  Input 


Depth  (fsw) : 

- -  j 

2000 

Initial/Desired  Temperatures  (°  F)  : 

90/80 

Initial/Desired  Relative  Humidities 

(%):  90/50 

Number  of  Divers: 

6 

Internal  Chamber  Volume  (£t^) : 

3300 

Recirculation  Flowrate  (cfm) : 

100 

Mission  Duration  (hrs) : 

240 

Diver  Working  Conditions: 

Light  Work 

Therefore,  the  ideal  gas  mixture  was  determined  to  be 
0.57%  O2  and  99.43%  He. 

3)  The  minimum  and  maximum  oxygen  percentages  were 
calculated  by  equations  2.7  and  2.8. 


( - )  (100%) 

/  903.6  . 

14.696 


%0,  =  0.26% 


%Oj  =  ( - (100%) 

,  903.6  > 

'  14.696  ' 


%0,  =  0.98% 
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since  only  mixtures  within  these  levels  are  ^cc^■  '.o, 

air  could  not  be  used  and  a  0-5%  Oj  and  99. 
gas  mixture  was  selected.  Therefore,  '’a. 
scrubbers  were  used  and  a  computation  of  free  ve  ij  ilation 
was  not  needed. 

4)  The  pressure  and  depth  ranges  for  the  gas  mixture 
were  calculated  from  equations  2.9  -  2.12. 

(14.696) 

^’mln  =  47  0.3  psia 


=  (■^)  (14.696) 
=  176  3.5  psia 


=  [(^^>(33.066)  -  33.066] 


=  1025.0  fSW 


cUx  =  [(4?^)  (33.066)  -  33.066] 

"  3934.9  fsw 

5)  Since  a  helium-oxygen  mixture  was  used  with  carbon 
dioxide  scrubbers,  the  amount  of  pure  oxygen  injected 
into  the  chamber  was  calculated  by  equation  2.22. 
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02c  “  0.0247  SCFM  for  light  work 


O2  =  (0.0247)  (6) 
O,  =  0.14  82  SCFM 


6)  The  amount  of  carbon  dioxide  produced  and  the 
scrubber  requirements  were  calculated  by  equations  2.17 
and  2.23  -  2.24. 


irt  =  (0.1137)  (0.9)  (0.1482) 
w  =  0.01517  Ib^/min 


=  (0.0157)  (240)  (60) 
iibco,  =  218.4  lb 


7fi  -  /  218.4\ 

-  ‘0:245’ 


lbs  =  891.4  lb 


7)  The  time  required  to  reach  95%  steady-state  was 
determined  from  equation  2.21  where  was  determined 


from  2.14. 


„  _  [  (0.01517)  (35.11)  (549.67)] 

~  10.58 

*  27  .7  ft^/min 


t  = 


(-0.0167)  ( 


3300.1  f  (0.2863) (27 .7)  , 

27 .7  ^  ^  (0.0152)  (549.67)  ^ 


t  =  5.96  hrs 


8)  Next,  the  diving-gas  properties  were  determined. 
Molecular  weight  was  calculated  from  equations  2 .25-2 .21 . 


0.5 

100 


0.005 


.  _  99.5 

'  100 


0.995 


MW  •  [  (0.005)  (32.00)  --  (0 . 99  5)  (4 . 0026  )  ] 
MW  ^  4.143  lb/ lb-mole 


Specific  heat  was  calculated  from  equation  2.29. 
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Cp  =  t  (0.005)  (6.95)  +  (0.995)  (4 .965)  ] 
Cp  =  4.97  5  BTU/ lb-mole 


Density  was  calculated  from  equations  2.35  -  2.38. 


B  =  ((11.44) (0.995)^  -  (8.8) (0.995) (0.005) 
-  (16.24)  (0.005)2] 

B  =  11.2187  cm^/g-mole 


B  =  (0.0160) (11.2817) 

B  =  0.1805  ft^/ lb-mole 


^  _  (1544)  (549.67) 

(903.6)  (144) 


K  =  6.5226 


_  6.5226  +  v/6. 52262  +  4  (6.5226)  (0.1805 
"  “  2 
=  0.6185  fc^/ lb-mole 


=  4-143 

^  6.6984 

p  =  0.6185  Ib/ft^ 
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9)  The  psychrometric  properties  were  calculated  from 
equations  2.39  -  2.48.  The  psychrometric  chart  created 
for  this  case  is  presented  as  Figure  4.2. 

Pg  =  0.69838  psia  (from  Table  2-7) 

Vj^  =  0.01610  cu  ft/lb  (from  Table  2-8) 

=  0.2880  cu  ft/mole 

.  _  [(0.69838)  (0.2880)  (903.6)  (144)] 

»  ■  [(1544)  (549.67)3 

APg  =  0.0308  psia 


Pg  =  0.69838  +  0.0308 
P^  =  0.7292  psia 


^  _  0.7297 _ 

”  903.6  -  0.7297 

X„g  =  0.00081  moles  water  vapor /mole  dry  gas 


(X^g  is  shown  as  point  A  on  Figure  4.2) 

=  (0.9)  (0.00081) 

=  0.00072  moles  water  vapor /mole  dry  gas 
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(X„i  is  shown  as  point  B  on  Figure  4.2) 


=  (0-5)  (0.00081) 

X^2  0.00041  moles  water  vapor/mole  dry  gas 


(X„2  is  shown  as  point  D  on  Figure  4-2) 

H  =  (20,000) (0.00081) 

H  =  16.2  BTU/mole  dry  gas 


Hi  =  (09) (16.2) 

Hi  =  14.6  BTU/mole  dry  gas 


Hj  =  (0.5) (16.2) 

Hj  =  8.1  BTU/mole  dry  gas 


10)  The  maximum  coil  temperature  was  determined  from  the 
psychrometric  chart  as  the  intersection  of 
100%  saturation  curve.  This  point  is  at  approximately 
70°  F  and  is  shown  as  point  C  on  Figure  4.2. 

Tj,  =  70°  F 


59 


11)  The  energy  requirements  were  calculated  from 
equations  2.51  -  2.56. 


-  (4.975)  (90)  ■*-  14.6] 

‘  ~  4.143 

=  111.6  BTU;ib 


_  [(4.975)  (70)  ■>-  8.1] 

4.143 


=  86 .0  BTU/lb 


=  (0.6185) (100) (111 .6  -86.0) 
=  1583 . 4  BTU/min 


Ah  = 


[4.975)  (80  -70)  ] 
4 . 143 


Ah  =  12.0  BTU/lb 


=  (0.6185'  (100)  (12.0) 
=  7  42.2  BTU/min 


Of.  =  1583 .4  +  742 . 2 
0^  =  2325.6  BTU/min 
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12)  The  moisture  removal  rate  was  calculated  from 
equations  2.57  -  2.59. 


=  (0.00072)  (  ) 

^  4.143 

=  0.00313  Ibjlb^ 


=  (0.00041)  (i§^^°^) 
^  4.143 

=  0.0017  8  Ibjlba 


MRR  =  (0.00313  -  0.00178)  (100)  (0.S185)  (1440) 
MRR  =  120.1  lb/ day 


The  actual  program  printout  is  provided  as  Table  4-2.  The 
slight  discrepancies  between  the  example  and  the  actual  output 
were  caused  by  round-off  error  obtained  while  working  rhe 
example. 


Table  4-2  Output  for  Hyperbaric  Chamber  Simulation 


HYPERBARIC  CHAMBER  SIMULATION 
Input  Conditions: 

Pressure  (psia) : 

Initial/Desired  Temperature  {°  F)  : 
Initial/Desired  Relative  Humidity  (%): 

#  of  Divers /Working  Conditions: 

Chamber  volume  (ft^) : 

Flow  Rate  (ft^/min): 

Mission  Duration  (hrs) : 

Output  Conditions: 

Diving-Gas  Mixture  (%He/%02)  : 

Oxygen  Required  (SCFM) : 

CO2  Produced  (lbs): 

Sodasorb  Required  (lbs): 

Time  to  95%  Steady-State  (hrs) : 

Molecular  Weight  ( lb/ lb-mole) : 

Gas  Mixture  Density  (Ib/ft^): 

Maximum  Coil  Temperature  (°  F)  : 

Dehumidif iration  Energy  (BTU/min) : 

Reheat  Energy  (BTU/min) ; 

Total  Energy  Requirement  (BTU/min) : 
Moisture  Removal  Rate  (Ib/day); 


903,6 

90.0/80.0 

00.0/50.0 

6.0/2 

3300.0 

100.0 

240.0 

99.5/0.5 

0.148 

218.7 

892.7 
6.0 

4.143 

0.6184 

70.2 
1563 . 5 

724 . 3 
2287 . 8 

125.2 


Exa^jDla  2: _ 8  ibi*arine  Diver  Lock-out  Chamber 


'^his  examp^d  sinulatas  a  diver  lock-out  chamber  aboard  a 
stibma'rine  used  to  decompress  divers  after  a  mission.  The 
input  conditions  for  this  example  are  given  in  Table  4-3. 


Table  4-3  Diver  Lock-out  Chamber  Simulation  Input 


Depth  (fsw) : 

400 

Initial/Desired  Temperatures  (°  F) 

85/78 

Initial/Desired  Relative  Humidities  (%) : 

90/60 

Number  of  Divers: 

2 

Internal  Chamber  Volume  (ft^) : 

100 

Recirculation  Flowrate  (cfm) : 

10 

Mission  Duration  (hrs) : 

4 

Diver  Working  Conditions: 

Resting 

The  calcvilation  steps  are  demonstrated  below: 

1)  Since  a  working  depth  was  input,  it  was  converted  to 
a  pressure.  By  equation  2.1, 


P  = 


,  (400  33.066) 

33.066 


) (14.696) 


P  =  192,5  psia 


2)  The  ideal  gas  mixture  for  this  depth  was  calculated 
by  equations  2.5  and  2.6. 
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%02 

%02 


( 


2 


0.35 
.  192.5 
14.696 

.67% 


) 


) (100%) 


%He  =  100  -  2.67 
%He  =  97 .33% 

► 

Therefore,  the  ideal  gas  mixture  was  determined  to  be 

2.67%  O,  and  97.33%  He. 

►  ^ 

3)  The  minimum  and  maximum  oxygen  percentages  were 
calculated  by  equations  2.7  and  2.8. 


%02 


•in 


%02 


■in 


, _ 0.16 

.  192.5 
14.696 

1.2% 


) 


) (100%) 


» 


I 


0  . 

192 . 5 
14.696 


(100%) 


4 . 58% 


i 


Since  only  mixtures  within  those  levels  air  acceptable, 
air  could  not  be  used  and  a  3%  O2  and  97%  He  mixture  was 
selected.  Therefore,  carbon  dioxide  scrubbers  were  used 
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■iMllllfli 


and  a  computation  of  free  ventilation  rate  was  not 
needed. 

4)  The  pressure  and  depth  ranges  for  the  gas  mixture 
were  calculated  from  equations  2.9  -  2.12. 


P.i„  '  (^><14.696) 
f.ln  -  '78-4 


(14.696) 
=  293.9  psia 


<4in  =  (33.066)  -  33.066] 

=  143.3  fsw 

cUx  =  ((  4°-^  >  ^33.066)  -  33.066] 
3.0 

--  6  28.2  fsw 


5)  Since  a  helium-oxygen  mixture  was  used  with  carbon 
dioxide  scrubbers,  the  amount  of  pure  oxygen  injected 
into  the  chamber  was  calculated  by  equation  2.22. 


0.5  SLM  or  0.0177  SCFM 


(0.0177) (2) 


O,  =  0.0354  SCFM 


6)  The  amount  of  carbon  dioxide  produced  and  the 
scrubber  requirements  were  calculated  by  equations  2.17 
and  2.23  -  2.24. 


m  =  (0.1137)  (0.9)  (0.0354) 
m  =  0.0036  Ib^/min 


Ib^  =  /0.0036)  (4)  (60) 
2b^^  =  0.86  Ib 


{  0 » ) 
'  0.245 


Ibj  -  3.5  Ib 


7)  The  time  required  to  reach  95%  steady-state  was 
determined  from  equation  2.21  where  was  determined 

from  2.14. 

V  =  [  (0.0036)  (35.11)  (544.67)  ] 

~  10.58 

=  6.51  f  t^/n\in 
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t 


(-0.0167)  ( 


100  nnfi-  (0-2863)  (6.51)  , 

6.51  (0.0036) (544.67) 


t  =  0.77  hrs 


8)  Next,  the  diving-gas  properties  were  determined. 
Molecular  weight  was  calculated  from  equations  2.25 
2.21 . 


3.0 

100 


0.03 


X  = 


100 


=  0.97 


[  (0.03)  (32.00)  ♦  (0.97)  (4 .0026)  ) 
MW  =  4.843  lb/ lb-mole 


Specific  heat  was  calculated  from  equation  2.29. 


Cp  =  [ (0.03)  (6.95)  +  (0.97)  (4.965)] 
Cp  =  5.025  BTU/ lb-mole 


Density  was  calculated  from  equations  2.35  -  2.38. 


B  =  [(11.44)  (0.97)2  -  (8.8)  (0.97)  (0.03) 
-  (16.24)  (0.03)2] 

B  =  10.49  32  cm^/g-mole 


B  =  (0.0160) (10.492) 

B  =  0.1679  ft^ lb-mole 


^  _  (1544)  (544.67) 

(192.5)  (144) 

K  =  30.34 


V  =  30.34  »  ^/30.542  4  (30.54)  (07l679) 

m  2 


V„  =  30.51  fty  lb-mole 


p  =  i-842 

^  30.51 

p  =  0.1587 

9)  The  psychrometric  properties  were  calculated  from 
equations  2.39  -  2.48.  The  psychrometric  chart  created 
for  this  case  is  presented  as  Figure  4.3. 

Pg  =  0.59604  psia  (from  Table  2-7) 

Vj^  =  0.01609  cu  ft/lb  (from  Table  2-8) 

=  0.2896  cu  ft/mole 
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Moisture  Mole  Fraction 


(moles  water  vapor/mole  dry  gas) 
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Div  FOiIb  Teiiipfrat  ure  (P) 

e  -l.d  l\s>chromet!'ir*  Chart  for  400  (fsw) 


^p 


3 


[(0,59604)  (0.2896)  (192.5)  (144)] 
(  (1544)  (544.67) ] 


LPg  ~  0.00569  psia 


P,  =  0.59604  +  0.00569 
Pg  =  0.60173  psia 


^  0.6C173 _ 

'  192.5  ••  0.60173 

=  0,00314  noles  water  vapor/ mole  dry  gas 


(X^^g  is  shown  as  point  A  on  Figure  4.3) 

X^^  =  (0.9)  (0.00314) 

^wi  ~  0<  00282  moles  water  vapor /mole  dry  gas 


is  shown  as  point  B  on  Figure  4.3) 


X^2  =  (0.6)  (0,00314) 

X^2  ~  0.00188  moles  water  vapor / mole  dry  gas 
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(X„2  is  shown  as  point  D  on  Figure  4.3) 


H  =  (20,000)  (0.00314) 

H  =  62.71  BTU/mole  dry  gas 


=  (0.9)  (62.71) 

=  56.44  BTU/mole  dry  gas 


H2  =  (0.6)  (62.71) 

=  37 .63  BTU/mole  dry  gas 


10)  The  maximum  coil  temperature  was  determined  from  the 

psychrometric  chart  as  the  intersection  of 

100%  saturation  curve.  This  point  is  at  approximately 

I 

70°  F  and  is  shown  as  point  C  on  Figure  4.3. 

T^  =  70°  F 

^  11)  The  energy  requirements  were  calculated  from 

equations  2.51  -  2.56. 

>  „  =  [(5.025)  (85)  56.44] 

"  4.843 

=  99 . 85  BTU/lb 
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I 


=  80.40  BTU/lb 

=  (0.1587) (10) (99.85  -  80.40) 
Od  =  30.9  BTU/min 


Ah  -  [  (5.025)  (78  -  70)  ) 

4.843 

Ah  =  8.3  BTUjlb 

0^  =  (0.1587) (10) (8.3) 
C?,  »  13.2  BTU/min 


=  30.9  +13.2 
C>t  *44.1  BTU/min 


12)  The  moisture  removal  rate  was  calculated  from 
equations  2.57  -  2.59. 


=  (0.00282)  ( 

^  4.843 

=  0.01049  Ibjlb^ 


=  (0.00188)  (  ) 
^  4.843 

=  0.00699  Ibjlb^ 


MRR  =  (0.01049  -  0.00699)  (10)  (0.1587)  (1440) 
MRR  =  8.0  lb/ day 


The  actual  program  printout  is  provided  as  Table  4-4.  The 
slight  discrepancies  between  the  example  and  the  actual  output 
were  caused  by  round-off  error  obtained  while  working  the 
example. 
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Table  4-4  Output  for  Diver  Lock-out  Chamber  Simulation 


DIVER  LOCK-OUT  CHAMBER  SIMULATION 


Input  Conditions; 

Pressure  (psia) : 

Initial/Desired  Temperature  (°  F) : 
Initial/Desired  Relative  Humidity  (%): 
#  of  Divers/Working  Conditions: 
Chamber  Volume  (ft^) : 

Flow  Rate  (ft^/min) : 

Mission  Duration  (hrs) : 

Output  Conditions: 

Diving-Gas  Moisture  (%He/%02): 

Oxygen  Required  (SCFM) : 

CO2  Produced  (lbs) : 

Sodasorb  Required  (lbs): 

Time  to  95%  Steady-State  (hrs) ; 
Molecular  Weight  ( lb/ lb-mole) : 

Gas  Mixture  Density  (Ib/ft^); 

Maximum  Coil  Temperature  (°  F)  ; 
Dehumidification  Energy  (BTU/min) : 
Reheat  Energy  (BTU/min) : 

Total  Energy  Requirement  (BTU/min) : 
Moisture  Removal  Rate  (Ib/day) ; 


192.5 

85.0/78.0 

90.0/60.0 

2.0/1 

100.0 

10.0 

4.0 

97.0/3.0 
0.035 
0.9 
3.5 
0.8 
4.843 
0.1587 
70.0 
30.9 
13.2 
44 . 1 
8 . 0 
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Example  3; _ Underwater  Habitat 


This  example  simulates  a  shallow  underwater  habitat.  The 
French  Conshelf  II  is  used  as  a  model.  Conshelf  II,  shown  as 
Figure  4.4,  was  a  shallow  water  habitat  operating  in  the  Red 
Sea  in  September  of  1963.  The  breathing  gas  was  air  and  the 
mission  duration  was  1  month.  Miller  and  Koblick  (1984) 
describe  this  mission  in  detail.  The  input  conditions  for 
this  example  are  given  in  Table  4-5.  The  calculation  steps 
are  demonstrated  below: 


Table  4-5  Underwater  Habitat  Simulation  Input 


Depth  (fsw) : 

50 

Initial/Desired  Temperatures  (°  F)  : 

95/80 

Initial/Desired  Relative  Humidities 

(%) :  95/50 

Number  of  Divers: 

4 

Internal  Chamber  Volume  (ft^): 

2000 

Diver  Working  Conditions: 

Moderate  Work 

1)  Since  a  working  depth  was  input,  it  was  converted  to 
a  pressure.  By  equation  2.1, 


P  = 


,  (50  »  33.066) 
33.066 


) (14.696) 


P  =  36.9  psia 


2)  The  Ideal  gas  mixture  for  this  depth  is  calculated  by 
equations  2.5  and  2.6. 


%C>j 


=  ( 


0.35 


%0,  = 


f  -3i.-  ?- ) 

14.696 

13.9% 


(100%) 


=  100  -13.9 
%He  =  86.1% 


3^ 


The  minitnum  and  maximum  oxygen  percentages  were 


calculated  by  equations  2.7  and  2.8. 


%02  =( - - (100%) 

(  36 .9  . 

'  14.696  ' 

%0,  =6.37% 


%o. 


-  /  -60 
'  /  36.9 
'  14.696 

=  23.9% 


) 


(100%) 


Since  any  mixture  within  these  levels  is  acceptable, 


air 


(21%  O2)  was  selected  us  the  diving-gas.  Therefore,  a 
free  ventilation  rate  was  used  to  keep  the  chamber  free 
of  carbon  dioxide  and  C02  absorbers  were  not  required. 


4)  The  free  ventilation  rate  at  depth  and  from  the 
surface  was  calculated  from  equations  2.14  -  2.17. 

^2c  ~  0.0530  SCFM  for  moderate  work 


CO,  =  (0.3)  (0.0530)  (4) 
CO,  =  0.1908  lb/ min 


*7?  =  (0.1137)  (0. 1908) 


ni  =  0.0217 


'air 


[  (0.0217)  (35.11)  (529  .67) ] 
10.58 


=  39.9  ft’/min 


S 


jir 


[(39.9)  (36.9)  (529.67)1 
[(14.696) (554.67) ) 


=95.8  SCFM 


5)  The  time  required  to  reach  a  95%  steady-state 
conditions  was  determined  from  equation  2.21. 


c 


(-0.0167) ( 


2000^1  Tj-  (0.2863)  (39.9) 
39.9  ’  (0.0217)  (554.67) 


] 


t  =  2 . 5  hrs 


6)  Next,  the  diving-gas  properties  were  determined. 
Since  the  diving-gas  was  air,  the  molecular  weights  and 
specific  were  29.0  Ib/lb-mole  and  6.95  BTU/lb -mole 

respectively.  Density  was  calculated  from  the  ideal 
equation  of  state,  equation  2.31. 


V8 


p 

p 


[(29.0)  (36.9)  (144) j 
[  (1544)  (554.67)1 

0.18002  Ib/ft^ 


7)  The  psychrometric  properties  were  calculated  from 
equations  2.39  -  2.48.  The  psychrometric  chart  created 
for  this  case  is  presented  as  Figure  4.5. 

Pg  =  0.81564  psia  (from  Table  2-7) 

Since  the  operating  pressure  was  less  than  100  psia,  no 
correction  was  required  for  the  saturation  vapor 
pressure. 

^  _  0.81564 


36.9  -  0.81564 


-  0.0226  moles  water  vapor /mole  dry  gas 


(X„g  is  shown  as  point  A  on  Figure  4.5) 

X^^  =  (0.95)  (0.0226) 

=  0.0215  moles  water  vapor /mole  dry  gas 


(X^i  is  shown  as  point  B  on  Figure  4.5) 
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Moisture  Mole  Fraction 


(moles  water  vapor/mole  dry 


r-l 

O 

O 


X^2  ^  0.0113  moles  wacer  vapoz/mole  dry  gas 

(X^2  is  shown  as  point  D  on  Figure  4.5) 

H  =  (20,000) (0.0226) 

H  =  452.0  BTU/mole  dry  gas 


=  (0.95)  (452.0) 

=  429.4  BTU/mole  dry  gas 


H2  =  (0.50)  (452.0) 

=  226.0  BTU/mole  dry  gas 


8)  The  maximum  coil  temperature  was  determined  from  the 
psychrometr ic  chart  as  the  intersection  of  Xw2  and  the 
100%  saturation  curve.  This  point  is  at  approximately 
73.7°  F  and  is  shown  as  point  C  on  Figure  4.5. 

T^,  =-  73.7°  F 

9)  The  energy  requirements  were  calculated  from 


equations  2.51  -  2.56. 


(  (6.95)  (95)  ^  429.4] 


n 


H. 


29.0 

^37.6  BTU/lb 


[  (6  .95)  (73.7)  226 . 0] 

29.0 

25.5  BTU/lb 


C»d  =  (0.18C02)  (39.9)  (37.6  -  25.5) 
=  86 . 9  BTU/min 


=  t6.95)  (80  -  73.7)] 

29.0 

A/i  *  1.51  BTU/lb 


Or  =  (0.18002)  (39.9)  (1.51) 
Or  =  10.8  BTU/min 


Oc  =  86.9  +  10.8 
Ot  =  97  .7  BTU/min 


10)  The  moisture  removal  rate  was  calculated  from 
equations  2.57  -  2.59. 


(0.18002) ( 18- 0158^ 
29 . 0 


W,  = 

=  0.01336  Ibjlba 


Wj  =  (0.18002)  (  ) 

^2  =  0.00702  Ibjlb^ 


MRR  =  (0.01336  -  0.00702)  (39.9)  (0.18002)  (1440) 
MRR  =65.5  lb/ day 


The  actual  program  printout  is  provided  as  Table  4-6.  The 
slight  discrepancies  between  the  example  and  the  actual  output 
were  caused  by  round-off  error  obtained  while  working  the 
example. 


Table  4-6  Output  for  Underwater  Habitat 


UNDERWATER  HABITAT  SIMULATION 

"  1 

Input  Conditions: 

Pressure  (psia) : 

36.9 

Initial/Desired  Temperature  (°  F) : 

95.0/80.0 

Initial/Desired  Relative  Humidity  (%) : 

95.0/50.0 

#  of  Divers/Working  Conditions: 

4.0/3 

Chamber  Volume  (ft^): 

2000.0 

Output  Conditions: 

Diving-Gas  Mixture: 

AIR 

Free  Ventilation  Rate  at  Depth  (cfm) : 

39.9 

Free  Ventilation  Rate  (SCFM) : 

95.8 

Time  to  95%  Steady-State  (hrs) : 

1 

2.5 

Molecular  Weight  ( lb/ lb-mole) : 

29.0 

Gas  Mixture  Density  (Ib/ft^): 

0.18002 

Maximum  Coil  Temperature  (°  F) : 

73.8 

Dehumidification  Energy  (BTU/roin) : 

87.0 

Reheat  Energy  (BTU/min) : 

10.8 

Total  Energy  Requirement  (BTU/min) : 

97.8 

Moisture  Removal  Rate  (Ib/day): 

1 

65.4 

CHAPTER  V 


SDMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A  FORTRAN  computer  program  has  been  developed  to  simulate 
a  hyperbaric  chamber  environment.  Although  hyperbaric 
chambers  have  been  used  for  many  years,  increasing  complexity 
requires  innovative  design  procedures.  This  program  allows 
the  user  to  design  life  support  systems,  temperature  and 
humidity  control  equipment,  and  determine  energy  requirements 
for  hyperbaric  chambers,  submarine  diver  lockout  chambers,  and 
underwater  habitats.  By  using  this  program,  the  design 
engineer  can  simulate  a  variety  of  scenarios  and  determine  the 
design  criteria  for  a  specific  chamber. 

Conclusions t 

1)  Acceptable  diving-gas  mixtures  are  obtained  by  ensuring 
the  oxygen  partial  pressure  of  the  mixture  remains  between 
0.16  and  1.6  atm  for  short  periods  of  time  and  below  0.6  atm 
for  extended  operations. 

2)  In  order  to  eliminate  the  narcosis  effects  of  nitrogen 
on  the  diver,  helium-oxygen  mixtures  are  used  below  150  fsw. 

3)  The  partial  pressure  of  carbon  dioxide  (CO2)  must  be 
maintained  below  0.005  atm  during  extended  missions  in  order 
to  eliminate  discomfort  and  possible  injury  to  the  divers. 
This  can  be  accomplished  by  ventilation  of  the  chamber  or  by 
using  a  CO2  absorption  system.  Ventilation  is  only 
recommended  when  air  is  used  as  the  diving-gas  due  to  the  cost 


of  helium. 


4)  A  CO2  absorption  system  consists  of  canisters  filled 
with  sodasorb.  The  diving-gas  mixture  flows  through  these 
canisters  and  the  CO2  is  removed.  Although  this  program  can 
be  used  to  estimate  the  duration  of  the  CO2  canisters  if  the 
amount  of  sodasorb  contained  is  known,  a  chemical  detection 
system  should  always  be  used  to  determine  when  to  replace  the 
absorbent  material. 

5)  The  molecular  weight  and  specific  heat  at  constant 
pressure  of  the  diving-gas  mixture  can  be  determined  from  the 
component  mole  fractions. 

6)  The  density  of  helium-oxygen  mixtures  is  determined  from 
the  virial  equation  of  state.  Three  virial  coefficients  are 
required  to  accurately  calculate  density.  Although,  the  first 
two  virial  coefficients  have  been  accurately  determined  for 
helium  and  oxygen,  only  one  second  virial  coefficient  for 
helium-oxygen  mixtures  has  been  determined.  Additionally,  no 
third  virial  coefficients  for  helium-oxygen  mixtures  have  been 
determined  to  date.  Therefore,  extremely  accurate  density 
calculations  for  helium-oxygen  mixtures  are  not  obtainable  at 
this  time. 

7)  Psychrometric  charts  have  been  used  for  years  to 
determine  dehumidification  and  cooling  requirements. 
Relationships  were  developed  to  compute  moisture  mole 
fractions  allowing  a  computer  simulation  without  using 
the  psychrometric  charts.  This  program  is  easier  to  use  and 
more  accurate  than  psychrometric  charts. 


8)  Enthalpy  of  moisture  content  (H)  is  ;aasily  calculated 
once  X„g  is  determined.  Enthalpy  is  an  important 
thermodynamic  property  for  calculating  energy  requirements. 

9)  Dehumidification  and  energy  requirements  can  be 
determined  for  varying  conditions  if  the  initial  and  desired 
temperatures  and  relative  humidities  are  known.  However,  the 
total  energy  requirement  neglects  heat  loss  due  to  heat 
transfer  between  the  chamber  and  the  environment. 
Additionally,  the  coil  temperature  required  to  perform  the 
desired  dehumidification  and  the  moisture  removal  rate  can  be 
determined. 

RaooDvuandatlons : 

1)  Pressure  vessel  design  of  a  hyperbaric  chamber  could  be 
ivicorporated  into  a  computer  program.  NAVFAC  DM-39  (1982) 
provides  guidance  for  pressure  vessel  design  of  hyperbaric 
chambers  and  the  Pressure  Vessel  Handbook  (1977)  simplifies 
the  calculations  of  pressure  vessel  design  for  desired 
conditions. 

2)  Fire  protection  and  extinguishing  systems  along  with 
other  safety  design  criteria  could  be  incorporated  into  a 
computer  program.  Additionally,  electrical,  utility,  and 
communication  systems  could  be  considered.  NAVFAC  DM-39 
(1982)  provides  guidance  for  these  systems. 

3)  Different  diving-gas  mixtures  could  be  investigated  with 
additional  facility  requirements  incorporated  into  the 
program.  The  NOAA  Diving  Manual  (1991)  lists  some 


k 

alternatives  to  he  .'v^m  and  research  into  the  effects  and 
advantages  of  these  gases  is  ongoing. 

4)  In  order  to  accurately  calculate  the  density  of  diving- 
gas  mixtures,  the  second  and  third  virial  coefficients  need  to 
be  accurately  determined.  Although  these  coefficients  exist 
for  helium  and  oxygen,  they  have  not  been  accurately 
determined  for  the  mixtures.  Texas  A  &  M  University  has  the 

^  facilities  to  accurately  determine  these  coefficients  if 

funding  is  available. 

5)  Heat  transfer  between  the  hyperbaric  chamber  and  the 

^  environment  is  not  considered  in  this  program.  A  detailed 

study  of  materials  and  insulation  would  be  necessary  prior  to 
incorporating  it  into  a  computer  program. 

I 

I 

I 

I 

» 
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APPENDIX  A 


NOMENCLATURE 

B  -  second  virial  coefficient,  ft^/roole 
C  -  third  virial  coefficient,  ft^/mole^ 

C02p  -  carbon  dioxide  production,  SCFM 
Cp  -  specific  heat  at  constant  pressure,  BTU/ lb-mole 
D  -  fourth  virial  coefficient,  ft^/mole^ 
d  -  working  depth,  feet  of  seawater  (fsw) 

d^a^  -  maximum  allowable  depth  for  a  diving-gas  mixture,  fsw 
-  minimum  allowable  depth  for  a  diving-gas  mixture,  fsw 
H  -  molar  enthalpy  of  moisture  content  at  saturation, 
BTU/mole  dry  gas 

Ha  -  moist-gas-mixture  enthalpy  at  initial  temperature,  BTU/ib 
Hg  -  moist-gas-mixture  enthalpy  at  coil  temperature,  BTU/ lb 

-  molar  enthalpy  of  moisture  content  at  initial  conditions, 

BTU/mole  dry  gas 

H2  -  molar  enthalpy  of  moisture  content  at  desired  conditions, 
BTU/mole  dry  gas 

lb^o2  "  carbon  dioxide  produced,  lb 
Ibg  -  sodasorb  required  for  absorption,  lb 
m  -  mass  flowrate  of  carbon  dioxide,  Ib/min 
MRR  -  moisture  removal  rate,  lb/ day 

MW  -  molecular  weight  of  diving-gas  mixture,  lb/ lb-mole 
n  -  number  of  moles 

-  number  of  divers 
O2C  -  oxygen  consumption,  SCFM 

0 
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021-  -  oxygen  required  to  sustain  divers,  SCFM 
P  -  absolute  pressure,  psia 

^max  "  maximum  allowable  pressure  for  a  diving-gas  mixture, 
psia 

^min  ”  minimum  allowable  pressure  for  a  diving-gas  mixture, 
psia 

PPo2  “  partial  pressure  of  oxygen,  atm 

PCO2  -  partial  pressure  of  carbon  dioxide,  atm 

Pg  -  vapor  pressure  of  water  in  saturated  moist  air,  psia 

^^®ideai  “  percentage  of  helium  in  ideal  diving-gas  mixture 

^®2ideai  "  percentage  of  oxygen  in  ideal  diving-gas  mixture 

^®2max  “  maximum  acceptable  oxygen  percentage 

%02n,in  -  minimum  acceptable  oxygen  percentage 

-  energy  required  for  dehumidification,  BTU/min 
Qj.  -  energy  required  to  reheat  from  coil  to  desired 
temperature,  BTU/min 
Q^.  -  total  energy  required,  BTU/min 
R  -  gas  constant  for  CO2,  35.11  Ib/lbmole-R 
RH^j  -  desired  relative  humidity,  % 

RH^  -  initial  relative  humidity,  % 

Rg  -  respiratory  quotient 

Ry  -  universal  gas  constant,  1544  f t-lb/ Ibmole-R 

-  ventilation  rate  of  air  at  standard  conditions,  SCFM 
T  -  initial  temperature,  °  F 

t  -  time  to  reach  95%  steady-state  conditions,  hr 

Tj,  -  maximum  coil  temperature  for  dehumidification,  F 

Tj  -  desired  temperature,  °  F 
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t(j  -  mission  duration,  hr 
U  -  internal  energy,  BTU 
V  -  volume,  ft^ 

-  volumetric  flowrate  at  depth,  ft^/min 
Vj,j,  -  volume  of  the  chamber,  ft^ 

-  specific  volume,  ft^/lb 
~  molar  volume,  ft^/mole 

-  humidity  ratio  at  initial  relative  humidity,  Ib^^/lb^ 
W2  -  humidity  ratio  at  desired  relative  humidity,  Ib^^/lb,^ 
Xne  “  mole  fraction  of  helium 

Xq2  '  mo] e  fraction  of  oxygen 
X„g  -  mole  fraction  of  water  at  saturation, 
moles  water  vapor/mole  dry  gas 
-  moisture  mole  fraction  at  initial  conditions, 
moles  water  vapor/mole  dry  gas 
X^2  “  moisture  mole  fraction  at  desired  conditions, 
moles  water  A^apor/mole  dry  gas 


AFPOOIX  B 
SOURCE  CODE 
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MAJOR  REPORT  -  "A  COMPUTER  SIMULATION  OF  A  HYPERBARIC 

CHAMBER  ENVIRONMENT.” 

REAL  VAR,  DEPTH,  PRESS,  TFI,  TRI ,  TFD,  RHI,  RHD, 
+NUMDIV,  CHVOL,  FLOWRATE,  DUR,  VARl,  02C0NS,  PER02 , 
+PERHE,  02MIN,  02MAX,  DMIN,  DMAX,  PSIAMIN,  PSIAMAX,  VAIR, 
+SURAIR,  02REQ,  TIME,  MW,  CP,  DENS,  BC,  BF,  XWS,  XWl, 
+XW2,  H,  HI,  H2,  HA,  HC,  CT,  QD,  QR,  QT,  MRR,  AIR,  LBC02 , 
+LBS,  CASES 

VAR  =  A  VARIABLE  TO  DETERMINE  WHETHER  THE  USER  WISHES  TO 
CALCULATE  THE  GAS  MIXTURE  BASED  UPON  DEPTH  OR 
PRESSURE. 

DEPTH  =  A  WORKING  DEPTH  USED  TO  CALCULATE  THE  DIVING-GAS 
MIXTURE  IN  FEET. 

PRESS  =  PRESSURE  AT  DEPTH  IN  PSIA. 

TFI  =  THE  INITIAL  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 
FARENHEIT. 

TRI  =  THE  INITIAL  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 
RANKINE. 

TFD  =  THE  DESIRED  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 
FARENHEIT. 

RHI  =  THE  INITIAL  RELATIVE  HUMIDITY  IN  THE  CHAMBER. 

RHD  =  THE  DESIRED  RELATIVE  HUMIDITY  IN  THE  CHAMBER. 
NUMDIV  =  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

CHVOL  =  THE  INTERNAL  VOLUME  OF  THE  CHAMBER  IN  CUBIC 
FEET. 

FLOWRATE  =  THE  FLOWRATE  INSIDE  THE  CHAMBER  USED  FOR 

RECIRCULATING  THE  DIVING-GAS  MIXTURE  WHEN  C02 
SCRUBBERS  ARE  USED. 

DUR  =  THE  MISSION  DURATION  IN  HOURS. 

VARl  =  A  VARIABLE  USED  TO  DETERMINE  THE  WORKING 
CONDITIONS. 

02C0NS  =  OXYGEN  CONSUMPTION  IN  CUBIC  FEET  PER  MINUTE. 
PER02  =  THE  PERCENT  OF  OXYGEN  IN  THE  DIVING-GAS  MIXTURE. 
PERHE  =  THE  PERCENT  OF  HELIUM  IN  THE  DIVING-GAS  MIXTURE. 
02MIN  =  MINIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  TT.L 
DIVING-GAS. 

02MAX  =  MAXIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 
DIVING-GAS. 

DMIN  =  THE  MINIMUM  DEPTH  IN  FEET  FOR  WHICH  A  GIVEN 
DIVING-GAS  CAN  BE  USED. 

DMAX  =  THE  MAXIMUM  DEPTH  IN  FEET  FOR  WHICH  A  GIVEN 
DIVING-GAS  CAN  BE  USED. 

PSIAMIN  =  THE  MINIMUM  PRESSURE  IN  PSIA  FOR  WHICH  A  GIVEN 
DIVING-GAS  CAN  BE  USED. 

PSIAMAX  =  THE  MAXIMUM  PRESSURE  IN  PriA  FOR  WHICH  A  GIVEN 
DIVING-GAS  CAN  BE  USED. 

VAIR  =  VENTILATION  RATE  AT  DEPTH  IN  CUBIC  FEET  PER 
MINUTE. 

SURAIR  -  VENTILATION  RATE  AT  THE  SURFACE  IN  STANDARD 
CUBIC  FEET  PER  MINUTE. 
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02REQ  =  OXYGEN  REQUIRED  TO  SUSTAIN  DIVERS  IN  STANDARD 
CUBIC  FEET  PER  MINUTE  IF  C02  ABSORBERS  ARE  USED. 
TIME  =  THE  TIME  REQUIRED  FOR  THE  CHAMBER  TO  REACH  99% 
STEADY-STATE  CONDITIONS. 

MW  =  THE  MOLECULAR  WEIGHT  OF  THE  DIVING-GAS  MIXTURE 
(LB/ LB-MOLE) . 

CP  =  THE  SPECIFIC  HEAT  OF  THE  DIVING-GAS  MIXTURE 
(BTU/MOLE) . 

DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU  FT)  . 
BC  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  CM/G-MOLE) . 

BF  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  FT/LB-MOLE) . 

XWS  =  THE  MOISTURE  MOLE  FRACTION  OF  WATER  VAPOR  AT 
SATURATION  (MOLES  WATER  VAPOR/ MOLE  DRY  GAS)  . 

XWl  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  INITIAL  RELATIVE 
HUMIDITY. 

XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 
HUMIDITY. 

H  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  SATURATION 
(BTU/MOLE  DRY  GAS). 

HI  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 
INITIAL  R.H. 

H2  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 
DESIRED  R.H. 

HA  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  INITIAL 
CONDITIONS  (BTU/LB  DRY  GAS) . 

HC  =s  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  DESIRED 
CONDITIONS. 

CT  =  THE  COIL  TEMPERATURE  REQUIRED  FOR  DEHUMIDIFICATION 
IN  DEGREES  FARENHEIT. 

QD  =  THE  ENERGY  REQUIRED  FOR  DEHUMIDIFICATION  (BTU/MIN) . 
QR  =  THE  ENERGY  REQUIRED  FOR  REHEAT  (BTU/MIN). 

QT  =  THE  TOTAL  ENERGY  REQUIRED  FOR  DEHUMIDIFICATION  AND 
REHEAT  (BTU/MIN) . 

MRR  =  THE  REQUIRED  MOISTURE  REMOVAL  RATE  FOR  THE 
DEHUMIDIFICATION  PROCESS  (LB/DAY) . 

LBC02  =  AMOUNT  OF  C02  PRODUCED  (LBS) . 

AIR  =  A  VARIABLE  TO  DETERMINE  WHETHER  OR  NOT  AIR  IS 
BEING  USED  AS  THE  DIVING-GAS. 

LBS  =  AMOUNT  OF  SODASORB  REQUIRED  FOR  C02  ABSORPTION 
(LBS) . 

CASES  =  THE  NUMBER  OF  CASES  TO  BE  CALCULATED. 


THIS  PORTION  OF  THE  PROGRAM  WILL  ALLOW  THE  USER  TO  INPUT 
THE  NECESSARY  PARAMETERS  FROM  A  DATA  FILE. 

OPEN  (UNIT=3 , FILE='VAR.DAT' , STATUS= ' OLD ' ) 

OPEN  (UNIT=4 , FILE= 'OUTPUT. DAT' , STATUS= ' UNKNOWN ' ) 

READ  (3,*)  CASES 

DO  500  I  =  1, CASES 
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READ  (3,*)  VAR,  DEPTH,  TFI ,  TFD,  RHI,  RHD 
READ  (3,*)  NUMDIV,  CHVOL,  FLOWRATE,  DUR,  VARl 


IF  (VAR.EQ.O)  THEN 

CALL  CONVERT (DEPTH,  PRESS) 
GO  TO  10 
ELSE 

PRESS=DEPTH 
GO  TO  10 
ENDIF 

10  TRI=TFI+459.67 

IF  (VARl.EQ.l)  THEN 
O2CONS=0.5/28. 3 
ELSEIF  (VAR1.EQ.2)  THEN 
02CONS=0.7/28.3 
ELSEIF  (VAR1.EQ.3)  THEN 
02C0NS=1.5/28.3 
ELSEIF  (VAR1.EQ.4)  THEN 
O2CONS=2.0/28.3 
ELSE 

O2CONS=3.0/28.5 

ENDIF 


THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  IDEAL  AND 
ACCEPTABLE  DIVING-GAS  MIXTURES. 

CALL  GASMIX  ( PRESS, 02MIN,02MAX) 

IF  (02MAX.GE. 21.0)  THEN 
AIR=1 

PER02=21. 0 
ELSE 


AIR=0 

IF  (02MAX.GE. 15.0)  THEN 
PERO2=15.0 

ELSEIF  (02MAX.GE. 10.0)  THEN 
PER02=10.0 

ELSEIF  (02MAX.GE. 5.0)  THEN 
PER02  =  5. 0 

ELSEIF  (02MAX.GE.3,0)  THEN 
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PER02=3.0 

ELSEIF  (02MAX.GE. 1.0)  THEN 
PER02=1.0 

ELSEIF  (O2MAX.GE.0.5)  THEN 
PER02=0. 5 
ELSE 

PER02=0. 1 
ENDIF 

PERHE=100 . 0-PER02 

CALL  GASRANGE  ( PER02  ,  DMIN ,  DMAX ,  PSI AMIN ,  PSIAM.\X ) 
ENDIF 


THIS  PORTION  OF  THE  PROGRAM  WILL  CALCULATE  THE  REQUIRED 
VENTILATION  RATE  WITH  OR  WITHOUT  CARBON  DIOXIDE 
ABSORBERS . 

IF  (AIR.EQ.l)  THEN 

CALL  FREERATE (02C0NS , NUMDIV, TRI , PRESS , VAIR , SURAIR) 
ELSE 

CALL  ABSRATE (NUMDIV, 02REQ, 02C0NS , DUR , LBC02 , LBS ) 

ENDIF 

CALL  STSTATE ( CHVOL , 02C0NS , NUMDIV , TRI , TIME) 

THIS  PORTION  OF  THIS  PROGRAM  CALCULATES  THE  PROPERTIES 
OF  THE  DIVING-GAS  MIXTURE. 

CALL  M0LWT(PER02  MW, AIR) 

CALL  SPECHEAT(PER02,CP,AIR) 

CALL  DENSITY ( PER02 , TRI , PRESS , MW , DENS ,  DC , BF , AIR) 

THIS  PORTION  OF  THE  PROGRAM  DETERMINES  THE  PSYCHROMETRIC 
PROPERTIES. 

CALL  PSYCHART ( TFI , TRI , PRESS , XWS . H ) 

CALL  XWSH (RHI , RHD , XWS , XWl , XW2 , H , HI , H2 ) 

CALL  COILTEMP(XW2 , PRESS, CT) 

CALL  MGMENTH(TFI,CT,CP,H1,H2,MW,HA,HC) 


t 


97 


C  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  ENERGY 

C  REQUIREMENTS  FOR  DEHUMIDIFICATION  AND  REHEAT  ALONG  WITH 

C  MOISTURE  REMOVAL  REQUIREMENTS. 

CALL  ENERGY ( DENS , AIR , VAIR, FLOWRATE , HA , HC , TFD , CT , CP ,  MW 
+QD,QR,QT) 

CALL  MOISTREM ( XWl , XW2 , AIR , FLOWRATE , MW , VAIR , DENS ,  MRR) 

C  THIS  PORTION  OF  THE  PROGRAM  PRODUCES  THE  OUTPUT. 

IF  (AIR.EQ.O)  THEN 
WRITE  (4,12) 

12  FORMAT  (IX,  '******************************************* 

+***★***■>«★**★★**********★***★•  ) 

WRITE  (4,15)  I 

15  FORMAT  (35X,'CASE  ',12) 

WRITE  (4,20) 

20  FORMAT  ( IX,  ' ******************************************* 

+****************************' ) 

WRITE  (4,25) 

25  FORMAT  ( 2 5X, '"HYPERBARIC  CHAMBER  SIMULATION"') 

WRITE  (4,30) 

30  FORMAT  ( IX,  ' ******************************************* 
+****************************') 

WRITE  (4,35) 

35  FORMAT  (5X, 'INPUT  CONDITIONS;') 

WRITE  (4,40)  PRESS 

40  FORMAT  ( 5X ,' PRESSURE  (PSIA) : ' , 24X, 1F6 . 1) 

WRITE  (4,45)  TFI,TFD 

45  FORMAT  ( 5X ,' INITIAL/ DESIRED  TEMP  ( DEG  F)  :  ' , 9X , 1F4 . 1 ,  ' / ' , 
+1F4 .1) 

WRITE  (4,50)  RHI,RHD 

50  FORMAT  (  5X,  '  INITIAL/DESIRED  REL  HUM  (%):', lOX , 1F4 . 1 ,'/' , 
+1F4.1) 

WRITE  (4,55)  NUMDIV,VAR1 

55  FORMAT  (5X,'#  OF  DIVERS/WORKING  CONDITIONS 7X , 1F4 . 1 , 

+'/' ,1F2.0) 

WRITE  (4,60)  CHVOL 

60  FORMAT  ( 5X ,' CHAMBER  VOLUME  (CU  FT) ; ' , 17X , IFG . 1 ) 
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WRITE  (4,70)  FLOWRATE 

70  FORMAT  (5X,'FLOW  RATE  (CU  FT/MIN) : ' , 18X , 1F6 . 1 ) 

WRITE  (4,75)  DUR 

75  FORMAT  (5X, 'MISSION  DURATION  (HRS) : ' , 17X , 1F6 . 1 ) 

WRITE  (4,100) 

100  FORMAT  (IX,'******************************************* 

-i-*  *******************  •  ) 

WRITE  (4,105) 

105  FORMAT  (5X, 'OUTPUT: ' ) 

WRITE  (4,110)  PERHE,PER02 

110  FORMAT  (5X, 'DIVING-GAS  MIXTURE  ( %HE/ %02 ) : ' , 9X , 1F4 . 1 , 
+'/ ' ,1F4.1) 

WRITE  (4,115)  DMIN,DMAX 

115  FORMAT  (5X, 'DIVING-GAS  MIXTURE  RANGE  ( FSW) : ' , 7X , 1F6 . 1 , 
+' / ' , 1F6. 1) 

WRITE  (4,120)  02REQ 

120  FORMAT  (5X, 'OXYGEN  REQUIRED  ( SCFM) : ' , 18X , 1F8 . 4 ) 

WRITE  (4,125)  LBC02 

125  FORMAT  (5X,'C02  PRODUCED  (LB) : ' , 22X, 1F6 . 1 ) 

WRITE  (4,130)  LBS 

130  FORMAT  (5X, ' SODASORB  REQUIRED; ', 22X, 1F6. 1) 

WRITE  (4,135)  TIME 

135  FORMAT  (5X,*TIME  TO  95%  STEADY  STATE  (HRS) : ' , 9X, 1F6 . 1) 
WRITE  (4,140)  MW 

140  FORMAT  (5X, 'MOLECULAR  WEIGHT  ( LB/MOLE) :', 13X , 1F8 . 3 ) 
WRITE  (4,145)  DENS 

145  FORMAT  (5X,'GAS  MIXTURE  DENSITY  (LB/CU  FT) ; ' , 9X , IFIO . 5) 
WRITE  (4,150)  CT 

150  FORMAT  ( 5X ,' MAXIMUM  COIL  TEMPERATURE  ( DEG  F)  :  '  , 7X , 1F6 . 1 ) 
WRITE  (4,155)  QD 

155  FORMAT  (5X, ' DEHUMIDIFICATION  ENERGY  (BTU/MIN) : ' , 6X , 
+1F6. 1) 

WRITE  (4,160)  QR 

160  FOR14AT  ( 5X,  '  REHEAT  ENERGY  (BTU/MIN)  16X,  1F6 . 1) 

WRITE  (4,165)  QT 

165  FORMAT  (5X, 'ENERGY  REQUIREMENTS  (BTU/MIM) : ' , lOX, 1F6 . 1) 


99 


WPITE  (4,170)  MRR 

FORMAT  (5X, 'MOISTURE  REMOVAL  RATE  (LB/DAY) 9X, 1F6 . 1) 
WRITE  (4,175) 

FORMAT  (IX, 'A****************************************** 
+****************************') 

ELSE 

WRITE  (4,200) 

FORMAT  (IX, '******************************************* 

+  *A****************-*r*********») 

WRITE  (4,205) 

FORMAT  ( 2 5X, '"HYPERBARIC  CHAMBER  SIMULATION"') 

WRITE  (4,210) 

FORMAT  (IX, '******************************************* 

+  ***********V;****************«) 

WRITE  (4,215) 

FORMAT  (5X, 'INPUT  CONDITIONS:') 

WRITE  .^RESS 

FORMAT  ,  iRESSURE  (PSIA) : ' , 24X, 1F6 . 1) 

WRITE  (4,225)  TFI,TFD 

FORMAT  (5X,  ' INITIAL/DESIRED  TEMP  (DEG  F) : ' , 9X , 1F4 . 1 ,  ' /  ' 
+1F4 . 1) 

WRITE  (4,230)  RHI,RHD 

FORMAT  (5X,  ' INITIAL/DESIRED  REL  HUM  (%):', lOX , 1F4 . 1 ,'/' , 
+1F4.1) 

WRITE  (4,235)  NUMDIV,VAR1 

FORMAT  (5X,'#  OF  DIVERS/WORKING  CONDITIONS 7X , 1F4 . 1 , 
,1F2.0) 

WRITE  (4,240)  CHVOL 

FORMAT  (5X, 'CHAMBER  VOLUME  (CU  FT) : ' , 17X, 1F6 . 1) 

WRITE  (4,245) 

FORMAT  (IX,  ' ******************************************1 

+****************************  '  j 

WRITE  (4,250) 

FORMAT  (5 X, 'OUTPUT: ' ) 

WRITE  (4,255) 

format  (5X, 'DIVING-GAS  MIXTURE 2 3X ,' AIR ' ) 


WRITE  (4,260) 

260  FORMAT  ( 5X DIVING-GAS  MIXTURE  RANGE  (FSW):',10X, 

+  '0.0/60.0'  ) 

WRITE  (4,265)  VAIR 

265  FORMAT  (5X,'FREE  VENTILATION  RATE  AT  DEPTH  (CFM):',3X. 
+1F6. 1) 

WRITE  (4,270)  SURAIR 

270  FORMAT  (5X,'FREE  VENTILATION  RATE  (SCFM) : ' , IIX, 1F6. 1) 
WRITE  (4,275)  TIME 

275  FORMAT  (5X,'TIME  TO  95%  STEADY-STATE  (HRS) : ' , 9X , 1F6 . 1) 
WRITE  (4,280)  MW 

280  FORMAT  (5X, ' MOLECULAR  WEIGHT  (LB/MOLE) 15X, 1F4 . 1) 
WRITE  (4,285)  DENS 

285  FORMAT  (5X, 'DENSITY  OF  AIR  (LB/CU  FT) :  '  ,  14X , IFIO . 5 ) 
WRITE  (4,290)  CT 

290  FORMAT  ( 5X ,' MAXIMUM  COIL  TEMPERATURE  (DEG  F)  :  '  , 7X, 1F6 . 1) 
WRITE  (4,295)  QD 

295  FORMAT  ( 5X, ' DEHUMIDIFICATION  ENERGY  (BTU/MIN) ; ' , 6X, 
+1F6. 1) 

WRITE  (4,300)  QR 

300  FORMAT  (5X, ' REHEAT  ENERGY  ( BTU/MIN) ;', 16X , 1F6 . 1 ) 

WRITE  (4,305)  QT 

305  FORMAT  ( 5X, ' ENERGY  REQUIREMENTS  (BTU/MIN) lOX , 1F6 . 1 ) 
WRITE  (4,310)  MRR 

310  FORMAT  ( 5X ,' MOISTURE  REMOVAL  RATE  (LB/DAY)  9X , 1F6 . 1 ) 
WRITE  (4,315) 

315  FORMAT  ( IX, ' ******************************************* 

■{-•k*-ii-k*'k**********ifk-k**-k******>  ) 

E?..  IF 

500  CONTINUE 


ooon  1.  n  Cj  oo 


SUBROUTINE  CONVERT (DEPTH, PRESS) 

THIS  SUBROUTINE  CONVERTS  A  WORKING  DEPTH  IN  FEET  TO  A 
PRESSURE  IN  PSIA. 

REAL  DEPTH,  PRESS 

DEPTH  =  A  WORKING  DEPTH  IN  FEET. 

PRESS  =  A  CONVERTED  PRESSURE  IN  PSIA. 

PRESS=( (DEPTH+33. 06633) /33. 06633) *14.696 

END 


SUBROUTINE  GASMIX ( PRESS , 02MIN , 02MAX ) 

THIS  SUBROUTINE  CALCULATES  AN  ACCEPTABLE  DIVING  GAS 
MIXTURE  BASED  UPON  A  MINIMUM  OXYGEN  PARTIAL  PRESSURE  OF 
0.16  ATM  AND  A  MAXIMUM  OXYGEN  PARTIAL  PRESSURE  OF  0.6 
ATM. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


REAL  PRESS,  PATM,  OIMIN,  02MIN,  OIMAX,  02MAX 


PRESS  =  THE  CHAMBER  PRESSURE  IN  PSIA. 

PATM  =  THE  CHAMBER  PRESSURE  IN  ATM. 

OJMIN  =  MINIMUM  ACCEPTABLE  OXYGEN  FRACTION  IN  THE 
DIVING-GAS. 

02MIN  MINIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 
DIVING-GAS. 

CiMAX  =  MAXIMUM  ACCEPTABLE  OXYGEN  FRACTION  IN  THE 
DIVING-GAS. 

02MAX  =  MAXIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 
DIVING-GAS. 


PATM=PRESS/14.696 
01MIN=0. 16/PATM 
02MIN=01MIN*100.0 
01MAX=0.6/PATM 
02MAX=^01MAX*  100 . 0 


END 
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C 


SUBROUTINE  GASRANGE (PER02 , DMIN , UMAX, PSI AMIN , PSIAMAX) 

THIS  SUBROUTINE  CALCULATES  A  DEPTH  AND  PRESSURE  RANGE 
FOR  A  GIVEN  DIVING-GAS  MIXTURE. 

REAL  PER02,  PMIN,  DMIN,  PMAX,  DMAX,  PSIAMIN,  PSIAMAX 

PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 
MIXTURE. 

PMIN  =  THE  MINIMUM  PRESSURE  (ATM)  BASED  UPON  AN  OXYGEN 
PARTIAL  PRESSURE  OF  0.16  ATM. 

DMIN  =  THE  MINIMUM  DEPTH  (FT)  FOR  A  GIVEN  DIVING-GAS. 
PMAX  =  THE  MAXIMUM  PRESSURE  (ATM)  BASED  UPON  AN  OXYGEN 
PARTIAL  PRESSURE  OF  0.6  ATM. 

DMAX  =  THE  MAXIMUM  DEPTH  (FT)  FOR  A  GIVEN  DIVING-GAS. 
PSIAMIN  =  THE  MINIMUM  PRESSURE  (PSIA)  FOR  A  GIVEN 
DIVING-GAS . 

PSIAMAX  =  THE  MAXIMUM  PRESSURE  (PSIA)  FOR  A  GIVEN 
DIVING-GAS. 

PMIN=16.0/PER02 
DMIN=PMIN*33. 06633-33 .06633 
PMAX=60.0/PER02 
DMAX=PMAX*33 .06633-33 .06633 
PSIAMIN=14 . 696*PMIN 
PSIAMAX=14 . 696*PMAX 

END 


SUBROUTINE  FREERATE  (02CONS ,  NUMDIV ,  TRI ,  PRESS ,  VAIR ,  SURAIR) 

THIS  SUBROUTINE  CALCULATES  THE  REQUIRED  VENTILATION  RATE 
FOR  THE  CHAMBER  WITHOUT  CARBON  DIOXIDE  ABSORBERS. 

REAL  02C0NS,  NUMDIV,  TRI,  PRESS,  VAIR,  SURAIR,  C02PR0D, 
+MFLOW 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


02C0NS  =  OXYGEN  C0NSU14PTI0N  IN  CUBIC  FEET  PER  MINUTE. 
NUMDIV  =  THE  N'JMBER  OF  DIVERS  IN  THE  CHAMBER. 

TRI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 
RANKINE. 

PRESS  =  THE  CHAMBER  PRESSURE  IN  PSIA. 

VAIR  =  THF  VF,^:7ILu^TI0N  RATE  AT  DEPTH  (CU  FT/MIN). 
SURAIR  T.HE  VENTILATION  RATE  AT  THE  SURFACE  (SCFM)  . 
C02PF0D  THE  TOTAL  AMOUNT  OF  CARBON  DIOXIDE  PRODUCED 
(CU  FT/MIN) . 

MFLOW  =  THE  MASS  FLOW  RATE  OF  CARBON  DIOXIDE  (LB/MIN) . 
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C02PROD=0 . 9*02C0NS*NUMDIV 
MFLOW=0. 11372474*C02PR0D 
VAIR=3 . 31817426*MFLOW*TRI 

•  SURAIR=(VAIR*PRESS*529.67) / (14.696*TRI) 

END 


C 


C 

C 


SUBROUTINE  ABSRATE ( NUMDIV , 02REQ , 02C0NS , DUR , LBC02 , LBS ) 

THIS  SUBROUTINE  CALCULATES  THE  REQUIRED  VENTILATION  RATE 
OF  THE  CHAMBER  IF  CARBON  DIOXIDE  ABSORBERS  ARE  UTILIZED. 

REAL  NUMDIV,  02REQ,  02C0NS,  DUR,  L3C02 ,  MFLOW,  LBS, 
+C02PR0D 


NUMDIV  =  THE  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

02REQ  =  THE  AMOUNT  OF  PURE  OXYGEN  REt-UIRED  (SCFM)  . 
02C0NS  =  THE  OXYGEN  CONSUMPTION  (SCFM) . 

DUR  =  THE  MISSION  DURATION  IN  HOURS. 

LBC02  =  THE  AMOUNT  OF  C02  PRODUCED  (LB) . 

MFLOW  =  THE  MASS  FLOWRATE  OF  CARBON  DIOXIDE  (LB/MIN) . 
C02PR0D  =  CARBON  DIOXIDE  PRODUCTION  (CU  FT/MIN) . 

LBS  =  AMOUNT  OF  SODASORB  REQUIRED  FOR  C02  ABSORPTION 
(LBS) . 

O2REQ=NUMDIV*02C0NS 
CO2PROD=0 . 9*02C0NS*NUMDIV 
MFLOW=0. 11372474*C02PR0D 
LBCO2=MFLOW*60 . 0*DUR 
LBS=LBCO2/0.245 

END 


SUBROUTINE  STSTATE ( CHVOL, 02C0NS , NUMDIV , TRI , TIME) 

THIS  SUBROUTINE  CALCULATES  THE  TIME  REQUIRED  FOR  THE 
CHAMBER  TO  REACH  A  95%  STEADY-STATE  CONDITION. 

REAL  CHVOL,  02C0NS,  NUMDIV,  TRI,  TIME,  C02PR0D,  MFLOW, 
+VAIR 

CHVOL  =  THE  INTERIOR  CHAMBER  VOLUME  (CUBIC  FEET) . 
02C0NS  =  OXYGEN  CONSUMPTION  (CU  FT/MIN) . 

NUMDIV  =  THE  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

TRI  =  THE  INITIAL  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 
RANKINE. 

TIME  =  THE  TIME  REQUIRED  TO  REACH  99%  STEADY-STATE 
(HRS) . 
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CO2PR00  =  THE  CARBON  DIOXIDE  PRODUCTION  (CU  FT/MIN) . 
MFLOW  =  THE  MASS  FLOWRATE  OF  CARBON  DIOXIDE  (LB/MIN) . 
VAIR  =  THE  VENTILATION  RATE  AT  DEPTH  (CU  FT/MIN) . 

COi:PROD=0 . 9*02C0NS*NUMDIV 
MFLOW=0. 1137 2474 *C02PR0D 
VAIR=3.31817426*MFLOW*TRI 

TIME=(-0.0167*CHVOL/VAIR) *LOG ( 1 . 0- ( 0 . 2863 *VAIR) 

+  / (MFLOW-TRI) ) 

END 


SUBROUTINE  MOLWT (PER02 , MW, AIR) 

THIS  SUBROUTINE  CALCULATES  THE  MOLECULAR  WEIGHT  OF  THE 
DIVING-GAS  MIXTURE. 

REAL  PER02,  MW,  X02 ,  XHE,  AIR 

PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 
MIXTURE. 

MW  =  THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  (LB/ LB-MOLE) . 
X02  =  THE  MOLE  FRACTION  OF  OXYGEN  IN  THE  MIXTURE. 

XHE  =  THE  MOLE  FRACTION  OF  HELIUM  IN  THE  MIXTURE. 

AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

IF  (AIR.EQ.l)  THEN 

HW=29.0 

ELSE 

X02=PER02/100. 0 
XHE=1.0-XO2 

MW=(X02*32.00)+(XHE*4.0026) 

ENDIF 

END 


SUBROUTINE  SPECHEAT ( PER02 , CP, AIR) 

THIS  SUBROUTINE  CALCULATES  THE  SPECIFIC  HEAT  OF  THE 
DIVING-GAS  MIXTURE. 

REAL  PER02,  CP,  X02 ,  XHE,  AIR 

PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 
MIXTURE. 
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CP  =  THE  SPECIFIC  HEAT  OF  THE  DIVING-GAS  MIXTURE 
(BTU/MOLE) . 

X02  =  THE  MOLE  FRACTION  OF  OXYGEN  IN  THE  MIXTURE. 

XHE  =  THE  MOLE  FRACTION  OF  HELIUM  IN  THE  MIXTURE. 

AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

IF  (AIR.EQ.l)  THEN 


CP=6.95 

ELSE 

X02=PER02/ 100.0 
XHE=1.0-X02 

CP=(X02*6.95)+(XHE*4.965) 

ENDIF 

END 


SUBROUTINE  DENSITY ( PER02 , TRI , PRESS , MW , DENS , BC , BF , AIR) 

THIS  SUBROUTINE  CALCULATES  THE  DENSITY  OF  THE  DIVING-GAS 
MIXTURE. 

REALPER02,  TRI,  PRESS,  MW,  DENS,  BC,  BF,  CONST,  MOLVOL, 
+AIR 

PER02  «  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 
MIXTURE. 

TRI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 
RANKINE. 

PRESS=  THE  INITIAL  CHAMBER  PRESSURE  IN  PSIA. 

MW  =  THE  MOLECULAR  WEIGHT  OF  THE  DIVING-GAS  MIXTURE 
(LB/ LB-MOLE) . 

DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU 
FEET) . 

BC  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  CM/MOLE) . 

BF  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  FT/MOLE) . 

CONST  =  RT/P 

MOLVOL  =  MOLAR  VOLUME  (CU  FT/ MOLE) 

AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 
IF  (AIR.EQ.l)  THEN 

DENS=(MW*PRESS*144 .0) / ( 1544 . 0*TRI ) 


X02=PER02/100. 0 
XHC=1.0-X02 


i 
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BC=(11.44*XHE**2) -(8.8*XHE*X02) -(16.24*X02**2) 

BF=BC*3 . 53 15E-05*454 . 0 

CONST=( 1544 . 0*TRI) / (PRESS *144 . 0) 

M0LV0L=(C0NST+SQRT( (CONST**2) +(4 . 0*CONST*BF) ) ) /2. 0 
DENS=MW/M0LV0L 

ENDIF 

END 


SUBROUTINE  PSYCHART (TFI , TRI , PRESS , XWS , H) 

THIS  SUBROUTINE  CALCULATES  THE  KOISTURE  MOLE  FRACTION 
AND  ENTHALPY  OF  MOISTURE  UTILT?''!.  .  THE  PSYCHROMETRIC 
CHART. 

REAL  TFI,  TRI,  PRESS,  XWS,  H,  PS,  DPS,  VL 

TFI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 
FARENHEIT, 

TRI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 
RANKINE. 

PRESS  =  THE  INITIAL  CHAMBER  PRESSURE  IN  PSIA. 

XWS  =  THE  MOISTURE  MOLE  FRACTION  (MOLES  WATER  VAPOR/MOLE 
DRY  GAS) . 

H  =  THE  ENTHALPY  OF  MOISTURE  CONTENT  (BTU/MOLE  DRY  GAS)  . 
PS  =  THE  SATURATION  PRESSURE  OF  WATER  (PSIA) . 

DPS  =  THE  CHANGE  IN  SATURATION  PRESSURE  OF  WATER  FOR 
PRESSURES  ABOVE  100  PSIA. 

VL  =  THE  SPECIFIC  VOLUME  OF  THE  SATURATED  LIQUID  (CU 
FT/LB) . 

THIS  PORTION  OF  THE  SUBROUTINE  CALCULATES  THE  SATURATION 
PRESSURE  OF  WATER  GIVEN  THE  DRY  BULB  TEMPERATURE. 

PS=(1.9285E-11*TFI**5)  +  (9.8683E-10*TFI**4)  -c 
+(3.3193E-07*TFI**3) +(1.9916E-05*TFI**2) + 
+(1.0488E-03*Tri) +2 . 2080E-02 


THIS  PORTION  OF  THE  SUBROUTINE  CORRECTS  THE  SATURATION 
PRESSURE  FOR  PRESSURES  ABOVE  100  PSIA  USING  POYNTING'S 
EQUATION . 

IF  (PRESS. GT. 100. 0)  THEN 

IF  (TFI. GT. 32.0. AND. TFI. LE. 52.0)  THEN 
VL=0. 01602 

ELSEIF  (TFI. GT. 52.0. AND. TFI. LE. 60.0)  THEN 
VL=0. 01603 

ELSEIF  (TFI. GT. 60.0. AND. TFI. LE. 66.0)  THEN 
VL=0 .01604 
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ELSEIF  (TFI. GT. 66.0. AND. TFI.LE. 71.0)  THEN 
VL=0, 01605 

ELSEIF  (TFI. GT. 71.0. AND. TFI.LE. 76.0)  THEN 
VL=0. 01606 

ELSEIF  (TFI. GT. 76.0. AND. TFI.LE. 80.0)  THEN 
VL=0. 01607 

ELSEIF  (TFI. GT. 80.0. AND. TFI.LE. 84.0)  THEN 
VL=0. 01608 

ELSEIF  (TFI. GT. 84.0. AND. TFI.LE. 87.0)  THEN 
VL=0. 01609 

ELSEIF  (TFI. GT. 87.0. AND. TFI.LE. 91.0)  THEN 
VL=0. 01610 

ELSEIF  (TFI. GT. 91.0. AND. TFI.LE. 94.0)  THEN 
VL=0. 01611 

ELSEIF  (TFI. GT. 94.0. AND. TFI.LE. 97.0)  THEN 
VL=0. 01612 

ELSEIF  (TFI. GT. 97.0. AND. TFI.LE. 100.0)  THEN 
VL=^0. 01613 

ELSEIF  (TFI. GT. 100.0)  THEN 
VL=0. 01614 

ENDIF 

DPS=(1.67876*PS*VL*PRESS) /TRI 
PS=PS+DPS 

GO  TO  40 

ELSE 

GO  TO  40 
ENDIF 

40  XWS=PS/ (PRESS-PS) 

H=20000.0*XWS 

END 


SUBROUTINE  XWSH  (RHI ,  RHD ,  XWS ,  XWl ,  XVv'2  ,  H ,  H 1 ,  H2 ) 

THIS  SUBROUTINE  CALCULATES  THE  INITIAL  AND  DESIRED 
MOISTURE  MOLE  FRACTIONS  AND  MOLAR  ENTHALPIES. 

REAL  RHI,  RHD,  XWS,  XWl,  XW2,  H,  HI,  H2 ,  RHI,  RH2 

RHI  =  THE  INITIAL  RELATIVE  HUMIDITY  IN  THE  CHAMBER  IN 
PERCENT . 

RHD  =  THE  DESIRED  RELATIVE  HUMIDITY  IN  THE  CHAMBER  IN 
PERCENT. 

XWS  =  THE  MOISTURE  MOLE  FRACTION  FOR  WATER  VAPOR  AT 
SATURATION. 
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XWl  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  INITIAL  RELATIVE 
HUMIDITY. 

XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 
HUMIDITY. 

H  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT 
SATURATION. 

HI  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 
INITIAL  R.H. 

H2  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AN  THE 
DESIRED  R.H. 

RHl  =  THE  INITIAL  RELATIVE  HUMIDITY  IN  THE  CHAMBER 
(FRACTION) . 

RH2  =  THE  DESIRED  RELATIVE  HUMIDITY  IN  THE  CHAMBER 
(FRACTION) . 


RH1=RHI/1C0.0 

RH2=RHD/ 100.0 

XW1=RH1*XWS 

XW2=RH2*XWS 

H1=RH1*H 

H2=RH2*H 


END 


SUBROUTINE  COILTEMP ( XW2 , PRESS , CT) 

THIS  SUBROUTINE  CALCULATES  THE  COIL  TEMPERATURE  REQUIRED 
FOR  DEHUMIDIFICATION  TO  A  DESIRED  RELATIVE  HUMIDITY. 

REAL  XW2,  PRESS,  PS,  CT 


XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 
HUMIDITY. 

PRESS  =  THE  INITIAL  CHAMBER  PRESSURE  (PSIA) . 

PS  =  THE  SATURATION  PRESSURE  OF  WATER  (PSIA) . 

CT  =  THE  COIL  TEMPERATURE  REQUIRED  FOR  DEHUMIDIFICATION. 

PS=(XW2*PRESS) / (1.0+XW2) 

CT=(2.4135E02*PS**5) -(8.4538E02*PS**4)  + 

+(1. 1635E03*PS**3) -(8. 1910E02*PS**2) + 
+(3.5407E02*PS)+7.0106 

END 


109 


oooon  no  n  o o o o o o o o n o o o n o  non 


C 


SUBROUTINE  MGMENTH (TFI , CT, CP, HI , H2 , MW , HA, HC) 

THIS  SUBROUTINE  CALCULATES  THE  MOIST-GAS-MIXTURE 
ENTHALPIES  AT  THE  INITIAL  CONDITIONS  AND  COIL 
TEMPERATURE. 

REAL  TFI,  CT,  CP,  HI,  H2,  MW,  HA,  HC 

TFI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 
FARENHEIT. 

CT  =  THE  REQUIRED  COIL  TEMPERATURE  FOR  DEHUMIDIFICATION 

(F)  . 

CP  =  THE  SPECIFIC  HEAT  OF  THE  GAS  MIXTURE  (BTU/MOLE) . 
HI  =  THE  MOLAR  ENTHALPY  OF  MOISTbLE  CONTENT  AT  THE 
INITIAL  R.H. 

H2  =  THE  MOLAR  ENTHALPY  OF  MOIOTUPE  CONTENT  AT  THE 
DESIRED  R.H. 

MW  =  THE  MOLECULAR  WEIGH'’’  OF  THE  MIXTURE  (LB/LB-MOLE)  . 
HA  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  INITIAL 
CONDITIONS. 

HC  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  DESIRED 
CONDITIONS. 

HA=( (CP*TFI)+H1) /MW 
nC=( (CP*CT)+H2; /MW 

END 


SUBROUT  I N  E  E.N  V  (DENS, AIR,  VAIR,  FLOWK  ATE ,  HA ,  HC ,  TFD ,  CT , 

+CP,MW,QD,QR,QT) 

THIS  SUBROUTIf^F.  CALCULATES  THE  ENERGY  REQUIRED  FOR 
DEHUMIDIFICATJON  AND  REHEAT. 

REAL  DENS,  AIR  VAIR,  FLOWRATE,  HA,  HC,  TFD,  CT,  CP,  MW, 
+QD,  DH,  QR,  QT 

DENS  =  THE  DENSITY  Of  THE  DIVING-GAS  MIXTURE  (LB/CU  FT) . 
AIR  =  A  VARIABLE  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 
VAIR  =  THE  VENTILATiO.V  OR  CIRCUSATIO/;  RATE  AT  DEPTH  (CU 
FT/MIN) . 

FLOWRATE  =  THE  CIRCULATION  RATE  (CU  FT/MIN) . 

C  HA  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  INITIAL 

C  CONDITIONS. 

C  HC  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  DESIRED 

C  CONDITIONS. 

C  TFD  =  THE  DESIRED  CHAMBER  TEMPERATURE  AFTER  REHEAT  (F) . 
C  CT  -  THE  COIL  TEMPERATURE  REQUIRED  FOR  DEHUMIDIFICATION 
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C  CP  =  THE  SPECIFIC  HEAT  OF  THE  DIVING-GAS  MIXTURE 

C  (BTU/MOLE) . 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  (LB/ LB-MOLE) . 

C  QD  =  THE  ENERGY  REQUIRED  FOR  DEHUMIDIFICATION  (BTU/MIN) . 
C  DH  =  THE  CHANGE  IN  ENTHALPY  FROM  THE  COIL  TEMPERATURE  TO 

C  THE  DESIRED  CHAMBER  TEMPERATURE  (BTU/MIN) . 

C  QR  =  THE  ENERGY  REQUIRED  FOR  REHEAT  (BTU/MIN) . 

C  QT  =  THE  TOTAL  ENERGY  REQUIRED  DEHUMIDIFICATION  AND 

C  REHEAT  (BTU/MIN) . 

IF  (AIR.EQ.O)  THEN 
VAIR=FLOWRATE 

ELSE 

VAIR=VAIR 

ENDIF 

QD=DENS*VAIR* (HA-HC) 

DH=CP*(TFD-CT) /MW 

QR=DENS*VAIR*DH 

QT=QD+QR 

END 

C  - 

SUBROUTINE  MOISTREM ( XWl , XW2 , AIR , FLOWRATE , MW , VAIR , 
+DENS,MRR) 

THIS  SUBROUTINE  CALCULATES  THE  MOISTURE  REMOVAL  RATE 
REQUIRED  FOR  DEHUMIDIFICATION. 

REAL  Wl,  W2,  XWl,  XW2,  AIR,  FLOWRATE,  MW,  VAIR,  DENS, 
+MRR 

C  Wl  =  THE  HUMIDITY  RATIO  AT  THE  INITIAL  RELATIVE 

C  HUMIDITY. 

C  W2  =  THE  HUMIDITY  RATIO  AT  THE  DESIRED  REIJVTIVE 

C  HUMIDITY. 

C  XWl  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  INITIAL  RELATIVE 
C  HUMIDITY. 

C  XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 
C  HUMIDITY. 

C  AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

C  FLOWRATE  =  THE  CIRCULATION  RATE  (CU  FT/MIN) . 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  DIVING-GAS  MIXTURE 

C  (LB/ LB-MOLE) . 

C  VAIR  =  THE  VENTILATION  RATE  AT  DEPTH  (CU  FT/MIN) . 

C  DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU  FT) . 

C  MRR  ^  THE  REQUIRED  MOISTURE  REMOVAL  RATE  FOR  THE 

C  DEHUMIDIFICATION  PROCESS  (LB/DAY) . 


Ill 


IF  (AIR.EQ.O)  THEN 
VAIR=FLOWRATE 

ELSE 

VAIR=%'AIR 

ENDIF 

W1=XW1* (18 . 0158/MW) 

W2=XW2* ( 18 . 0158/MW) 
MRR=(W1-W2) *VAIR*DENS*1440. 0 

END 


C 


